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20  ABSTRACT  (Continue  on  nnM  aide  It  miMWr  md  Identity  by  black  number) 

Free  electron  lasers  based  on  stimulated  scattering  from  relativistic  electron  beams  show  potential  for 
a new  class  of  coherent  radiation  sources  ranging  from  millimeter  wavelengths  down  to  optical  wave- 
lengths and  beyond  which  are  characterized  by  continuous  tunability,  very  high  power  levels  and  com- 
petitive efficiencies.  Here  we  present  a comprehensive  fully-relativistic  treatment  of  stimulated  back- 
scattering  from  an  electron  beam  with  streaming  velocitwfl^eP  The  incident  pump  wive  is  taken  to  bel 
a static,  circularly  polarized,  periodic  magnetic  field  with  wavenumber^fc^'the  frequency  of  the  scatteredf 
radiation  is ■$I0cknH  1 The  analysis  is  linear  in  the  scattered  wave  fields,  but  contains  the  pump 
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field  to  all  orders.  The  physical  mechanism  of  the  backscattering  is  provided  by  the  v X B ponderjno- 
tive  forces  acting  on  the  beam  electrons  in  the  presence  of  both  the  pump  and  the  scattered  waves. 

This  force  produces  a longitudinal  charge  bunching  which  in  turn  gives  rise  to  source  currents  for  the 
scattered  waves, 

The  ankiysisis  steady  state  but  spatially  bounded.  In  addition  to  the  ponderomotive  force,  other 

effects  which  are  incorporated  are  the  self-consistent  collective  scalar  potential,  and  free-streaming  terms 
which  arise  from  boundary  conditions.  Various  regimes  of  the  interaction  may  be  distinguished  accord- 
ing to  which  of  these  effects  are  dominant;  viz.,  (1)  a low  gain,  single  particle  regime  appropriate  to 
tenuous  cold  beams,  in  which  the  boundary  terms  lead  to  spatial  transients;  (2)  a high-gain,  single- 
particle  regime  of  a tenuous  thermal  beam,  (stimulated  Compton  scattering);  (3)  a high-gain  cold  beam 
regime.  This  last  regime  is  subdivided  into;  (a)  a weak-pump  limit,  in  which  the  longitudinal  density 
wave  satisfies  normal-mode  dispersion  and  the  collective  scalai  potential  dominates  the  ponderomotive 
potential;  and  (b)  a strong-pump  limit,  in  which  the  ponderomotive  foice  dominates. 

From  a heuristic  argument  based  on  the  fact  that  the  saturation  mechr.nism  for  each  of  these 
processes  is  the  trapping  of  the  beam  electrons  in  the  potential  wells  of  the  longitudinal  wave,  we 
estimate  the  saturated  amplitude  of  the  backscattered  wave  and  the  efficiency  of  energy  extraction  from 
the  beam.  In  addition,  we  find  the  integrated  gain  over  the  system  length  for  the  spatially  transient 
case  (1),  and  the  exponentiation  length  in  the  spatially  homogeneous  cases  (2)  and  (3a, b). 

Finally,  we  discuss  the  present  status  of  experimental  work  both  with  linear  accelerators  and  with 
intense  relativistic  electron  beams. 

Agreement  between  experimental  results  and  the  predictions  of  the  linear  theory  in  regimes  1 and 
3a  is  good.  The  other  regimes  have  not  yet  been  subjected  to  experiment  and  no  definitive  experimental 
test  of  the  nonlinear  theory  has  been  made  in  any  regime.  Theoretical  efficiencies  are  largest  for  case 
3b,  and  we  present  sample  calculations  of  gigawatt  lasers  with  1 1%  efficiency  at  a wavelength  of  1700  /am 
and  0.25%  efficiency  at  0.58  ttm.  Possible  techniques  for  substantially  improving  on  these  efficiencies 
are  also  discussed. 
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FREE  ELECTRON  LASERS  AND  STIMULATED  SCATTERING 
FROM  RELATIVISTIC  ELECTRON  BEAMS 


I.  INTRODUCTION 

The  possibility  of  developing  lasers  in  which  the  active  medium  is  a stream  of  free  elec- 
trons has  recently  evoked  much  interest.  The  potential  advantages  are  numerous  and  include 
both  continuous  frequency  tunability  ( e.g . through  the  variation  of  electron  energy)  and  very 
high  power  operation,  since  no  damage  can  occur  to  the  lasing  medium  as  happens  with  liquid 
or  solid-state  lasers.  Furthermore,  it  is  hoped  that  the  high  efficiencies  which  characterize  free 
electron  generators  of  microwave  radiation  can  eventually  be  realized  with  free  electron  lasers. 

In  the  past  few  years,  the  first  operation  of  high  peak  power  free  electron  lasers  has  been 
reported.  They  have  generated  seven  kilowatts  in  the  infrared  regime  (Deacon  et  al.,  1977) 
and  megawatt  power  levels  in  the  infrared  (Granatstein  et  al.,  1977).  Furthermore,  recent 
theoretical  results  indicate  that  much  larger  powers  should  be  achievable  The  operative 
mechanism  in  these  experiments  was  stimulated  backscattering  of  a low  frequency  pump  wave 
from  a relativistic  electron  beam,  leading  to  amplification  of  coherent  radiation  at  a high  fre- 
quency corresponding  to  a double  doppler  shift  of  the  pump  frequency. 

The  low  frequency  pump  wave  may  take  the  form  of  a static  spatially  periodic  magnetic 
field  or  a propagating  electromagnetic  wave.  In  the  case  that  the  pump  wave  is  a static  periodic 
magnetic  field  of  period  /,  the  frequency  of  the  backscattered  wave  will  be  given  by 
(1  + v./c)y2v.(2ir/l),  which  for  a highly  relativistic  beam  becomes  — 4iry?c// where  v.  is  the 
axial  velocity  of  the  electrons  and  y:  - (1  - v?/c2)~l/2.  If  the  pump  field  is  an  electromagnetic 
pump  at  frequency  <o„  and  wavenumber  k„  - a >Jc,  propagating  antiparallel  to  the  electrons,  the 
upshifted  backscattered  frequency  will  be  (1  + v./c)2y}o>0  = 4y.2a>„.  For  either  type  of 
Manuscript  submitted  Nov.  14,  1978. 
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incident  pump  wave  the  scattered  wave  frequency  can  be  varied  simply  by  changing  the  electron 
beam  energy  or  the  incident  pump  wavelength. 


Stimulated  scattering  of  photons  by  an  electron  ensemble  was  first  predicted  by  Kapitza 
and  Dirac  (1933).  The  first  analyses  of  high  frequency  radiation  from  a relativistic  electron 
beam  passing  through  a periodic  magnetic  field  (undulator)  was  made  by  Motz  (1951).  The 
proposal  for  generating  short  wavelength  radiation  by  stimulated  scattering  of  an  electromag- 
netic pump  wave  from  a relativistic  electron  beam  is  due  to  Pantell  et  al.  (1968). 


Experiments  demonstrating  stimulated  scattering  with  an  output  in  the  infrared  have  been 
carried  out  at  Stanford  University.  These  experiments  utilize  the  relatively  low  current  electron 
beam  in  a linear  accelerator  (/  = l A,  24  MV  < V < 43  MV).  They  are  in  the  aingle-particle 
scattering  regime  and  hence  do  not  involve  collective  electron  oscillations.  The  single  pass 
gain  is  low  and  the  process  is  strongly  dependent  on  the  finite  length  of  the  interaction  region. 
The  first  Stanford  experiment  (Elias  et  al.,  1976)  demonstrated  wave  amplification  at  a 
wavelength  of  10.6  jim  using  a static  periodic  2.4  kG  helical  magnetic  pump  with  a 3.2  cm 
period;  a 7%  increase  in  power  was  achieved  in  a 5.2  m interaction  length.  The  second  Stanford 
experiment  (Deacon  et  al.,  1977)  produced  laser  oscillations  in  an  optical  cavity  with  peak 
power  of  7 kW  at  a wavelength  of  3.4  fim;  0.01%  of  the  electron  beam  energy  was  converted 
into  radiation. 

The  analysis  by  Motz  (1951)  is  applicable  to  the  single-particle  scattering  process  in  the 
Stanford  experiments.  A subsequent  analysis  by  Sukhatme  and  Wolff  (1973)  included  the 
important  effect  of  the  finite  length  of  the  interaction  region.  The  treatment  of  Sukhatme  and 
Wolff  (1973)  was  quantum  mechanical,  as  were  analyses  by  Madey  (1971)  and  Madey  et  al. 
(1973).  However,  recent  analyses  (Hopf  et  al.,  1976a;  Kroll  and  McMullin,  1978)  have  shown 
that  a classical  treatment  of  the  stimulated  scattering  process  is  appropriate.  Additional, 
detailed  analyses  of  the  Stanford  experiments  have  also  appeared  in  the  literature  (Colson, 
1976,  1977;  Bernstein  and  Hirshfield,  1978;  Sprangle  and  Granatstein,  1978),  as  have  con- 
siderations of  extending  operation  to  shorter  wavelengths  up  to  the  x-ray  regime  (Hasegawa  et 
al.,  1976;  Elias  et  al.,  1976). 
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With  intense  relativistic  electron  beams,  collective  effects  play  an  important  role  in  the 
scattering  process.  The  theoretical  analyses  (Sprangle  and  Granatstein,  1974;  Sprangle  et  at., 
1975)  show  that  for  intense  magnetized  relativistic  electron  beams,  the  pump  wave  may  decay 
into  either  a space  charge  beam  mode  or  a cyclotron  beam  mode  depending  on,  among  other 
things,  the  pump  frequency.  The  axial  magnetic  field  can  be  used  to  increase  the  effective 
strength  of  the  pump  field  by  adjusting  the  electron  cyclotron  frequency  to  equal  roughly  the 
pump  frequency  in  the  beam  frame.  The  resonance  produced  in  this  way  can  substantially 
enhance  the  gain  of  the  scattered  wave. 

Experimental  and  theoretical  research  demonstrating  stimulated  scattering  in  the  submil- 
limeter regime  have  been  carried  out  at  the  Naval  Research  Laboratory.  In  experiments  using 
intense  relativistic  electron  beams  (/  — 30  kA,  1 MV  < V < 3 MV),  strong  submillimeter 
radiation  was  first  reported  by  Granatstein  et  at.  (1974).  Analyses  of  this  process  by  Sprangle 
and  Granatstein  (1974)  and  Sprangle  et  at.  (1975)  show  that  the  submillimeter  radiation  was 
generated  by  stimulated  Raman-type  scattering.  That  is,  one  of  the  scattered  waves  was  a beam 
space  charge  wave  enhanced  by  a resonance  between  the  pump  frequency  and  the  electron 
cyclotron  frequency.  A followup  experiment  designed  to  enhance  the  stimulated  Raman 
scattering  process  generated  a megawatt  of  radiation  in  a superradiant  oscillator  operating  at  a 
wavelength  of  400  /an  (Granatstein  et  at.,  1977). 

Experiments  at  Columbia  University  using  a static,  rippled,  magnetic-field  pump  with  out- 
puts in  the  centimeter  wave  regime  (Efthimion  and  Schlesinger,  1977)  and  in  the  millimeter 
wave  regime  (Marshal  et  at.,  1977;  Gilgenbach  et  at.,  1978)  clearly  distinguished  between  the 
cases  where  the  idler  was  a space-charge  mode  and  where  the  idler  was  a beam  cyclotron  mode. 
Efficiencies  as  large  as  0.2%  were  measured  in  these  experiments.  Recently  a joint  Columbia- 
NRL  experiment  (McDermott  et  at.,  1978)  has  reported  operation  of  a stimulated  Raman 
scattering  laser  which  employs  a quasi-optical  cavity;  laser  output  was  1 MW  at  400  jxm  and  line 
narrowing  to  Aw/w  - 2%  was  observed  (compared  with  Aw/w  >10%  in  the  superradiant  oscil- 
lator case) . 
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Theoretical  treatments  of  stimulated  Raman  scattering  in  plasma  were  first  applied  to 
laser-pellet  fusion  studies  as  well  as  to  ionospheric  experiments  (Silin,  1965,  1967;  DuBois  and 
Goldman,  1965,  1967;  Nishikawa,  1968;  Kaw  and  Dawson,  1971;  Perkins  and  Raw,  1971; 
Forslund  et  al. , 1973;  Manheimer  and  Ott,  1974;  Drake  el  al.,  1974).  Relativistic  analyses  with 
application  to  the  generation  of  short  wavelength  radiation  using  electron  beams  are  now  in  the 
literature,  treating  both  the  linear  regime  (Sprangle  and  Granatstein,  1974;  Miroshnichenko, 
1975;  Sprangle  et  al.,  1975;  Kroll  and  McMullin,  1978)  and  the  non-linear  regime  (Kwan  et  al., 
1977;  Sprangle  and  Drobot,  1978).  Efficiencies  greater  than  10%  have  been  calculated,  but  no 
efficiency  approaching  this  level  has  been  measured  to  date;  experiments  have  not  been  driven 
into  saturation.  However,  e-folding  lengths  on  the  order  of  ten  centimeters  have  been  calcu- 
lated and  correspond  closely  to  experimental  measurements  (Granatstein  et  al.,  1977). 


The  present  paper  confines  itself  to  free  electron  lasers  based  on  stimulated  scattering 
from  relativistic  electron  beams.  Other  mechanisms  involving  relativistic  electron  beams  which 
may  lead  to  free  electron  lasers  have  also  been  studied.  While  none  of  these  has  to  date  pro- 
duced the  high  power  in  the  submillimeter  and  infrared  regimes  that  characterize  the  stimulated 
scattering  work,  the  research  is  in  much  too  early  a stage  to  ignore  the  other  possible  mechan- 
isms. For  reference  we  cite  the  following  mechanisms  which  have  been  considered  as  the  basis 
for  free  electron  sources  of  coherent  high  frequency  radiation:  (1)  interaction  of  a relativistic 
electron  beam  with  periodic  slow  wave  structures  or  with  dielectric  media  (Coleman,  1961; 
Walsh  et  al.,  1977;  Schneider  and  Spitzer,  1974;  Gover  and  Yariv,  1978);  and  (2)  scattering  of 
electromagnetic  waves  from  the  front  of  a relativistic  electron  beam  or  from  a moving  ioniza- 
tion front  (Granatstein  et  al.,  1976;  Pasour  et  al.,  1977;  Buzzi  et  al.,  1977;  Lampe  et  al.,  1978). 


The  linear  regime  has  been  discussed  in  a survey  by  Granatstein  and  Sprangle  (1977)  and 
in  comprehensive  treatments  by  Kroll  and  McMullin  (1977)  and  Hasegawa  (1977).  The 
present  paper  presents  an  original  reformulation  of  stimuated  scattering  theory  that  attempts  to 
encompass  all  previous  treatments,  and  compares  the  theory  with  available  experimental  data. 
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In  addition,  we  present  a semi-qualitative  discussion  of  non-linear  saturation.  Also,  the  linear 
treatment  differs  from  previous  analyses  by  performing  a perturbation  expansion  about  an  exact 
electron  beam  equilibrium  in  a general  way. 

The  details  of  the  physical  mechanism  associated  with  FELs  depend,  to  some  extent,  on 
the  type  of  scattering  process  under  consideration,  which  in  turn  is  a function  of  the  system 
parameters,  such  as  beam  energy,  temperature  and  density,  the  length  of  the  interaction  region, 
strength  of  pump  field,  etc.  There  are  characteristics  of  the  physical  mechanism,  however, 
which  are  common  to  all  the  various  scattering  processes.  These  common  characteristics  can 
best  be  described  classically  in  the  beam  frame  of  reference. 

In  the  following  discussion,  quantities  in  the  beam  frame  will  be  written  with  primes.  In 
the  beam  frame  we  stipulate  that  the  existing  electron  equilibrium  is  perturbed  by  a low  fre- 
quency density  wave;  the  existence  of  this  wave  will  be  justified  later.  Only  waves  propagating 
along  the  z axis,  i.e.,  the  direction  of  the  beam  velocity  in  the  laboratory  frame,  will  be  con- 
sidered. This  longitudinal  perturbation  has  a frequency  and  wavenumber  given  by  (at',  k1). 
The  introduction  of  a large-amplitude  high-frequency  incident  electromagnetic  pump,  £',  at 
(a>o,  k'„)  forces  the  electrons  to  oscillate  at  a frequency  a >'  > > &»'  in  the  direction  along  £'  with 
a maximum  velocity  given  by  = |e|£'/(m„o>').  This  transverse  oscillation  velocity  v'Q1, 
perpendicular  to  k'„,  couples  to  the  density  wave,  thus  inducing  a transverse  current  at  fre- 
quency o/+  - <*>'  + t u'„  and  wave  number  k'+  = k1  + k'„.  This  current  now  generates  an  elec- 
tromagnetic wave  at  (a>‘+,  k'J.  The  generated  or  scattered  electromagnetic  field  consists  of 
backscattered  waves  propagating  antiparallel  to  the  incident  pump  wave.  Forward  scattered 
waves  are  also  induced,  but  will  not  be  considered  because  they  are  downshifted  in  frequency 
when  transformed  back  to  the  laboratory  frame.  The  pump  and  backscattered  wave  couple 
through  the  r'  x B'  term  in  the  Lorentz  force  equation,  resulting  in  a longitudinal  force  at 
(at',  k').  This  induced  longitudinal  force,  also  called  the  ponderomotive  or  radiation  pressure 


force,  will  reinforce  the  original  longitudinal  wave.  The  backscattered  electromagnetic  wave  is, 
therefore,  unstable  resulting  in  stimulated  emission  of  radiation. 


The  ponderomotive  potential  associated  with  the  (w , k)  wave,  may  be  comparable  to  or 
may  dominate  the  self  space  charge  potential  wave  set  up  by  the  beam.  If  the  frequency  of  the 
pump  wave  is  far  above  any  of  the  characteristic  frequencies  of  the  beam  we  find  that 
<a'  - - ck'0  and  o>+  - ck'+.  The  backscattered  wave  frequency  in  the  laboratory  is 
<o+  - y;(<u'  + v.k'J  — (1  + v./c)y.<u£  — (1  + v./c)2y2 <o„.  For  a highly  relativistic  beam  the 
backscattered  frequency  is  w+  ~ 4y.2a>„. 

In  Section  II  we  discuss  in  detail  the  physical  model  to  be  analyzed,  which  consists  of  a 
relativistic  electron  beam  of  arbitrary  intensity  entering  and  propagating  through  a static, 
periodic,  circularly  polarized  magnetic  pump  field.  The  general  analytic  formalism  is  presented. 
The  calculation  is  relativistically  covariant  and  consists  of  a Vlasov-Maxwell  perturbation 

| 

analysis.  The  perturbation  is  performed  about  an  exact  electron  beam  equilibrium  distribution, 
formed  from  constants  of  the  motion  in  the  magnetic  pump  field.  The  general  analysis  cul- 
minates in  expressions  for  the  driving  current  densities  which,  through  the  wave  equations,  are 
sources  of  the  backscattered  fields.  The  time  dependences  of  the  scattered  fields  are  harmonic, 
i.e.,  temporal  steady  state  conditions  exist.  Spatial  transient  effects  due  to  the  boundary  of  the 
interaction  region  at  z = 0 are  retained  by  choosing  the  scattered  fields  to  have  an  arbitrary  spa- 
tial dependence.  In  Section  III  we  obtain  spatial  growth  rates  and  estimates  for  the  saturation 
efficiencies  and  field  amplitudes  of  various  scattering  processes:  Case  (1)  is  a transient,  low- 
gain,  single-particle  regime  and  applies  to  the  recent  FEL  experiments  at  Stanford  University; 
Case  (2)  is  a high-gain,  single-particle,  thermal-beam  regime,  which  is  equivalent  to  the  well 
known  stimulated  Compton  scattering  process;  Case  (3a)  is  a high-gain,  cold  beam,  weak  pump 
scattering  regime  where  collective  space  charge  effects  play  a dominant  role;  Case  (3b)  is  the 
strong  pump  limit  of  Case  (3a).  In  Case  (3b)  collective  space  charge  effects  are  dominated  by 
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the  ponderomotive  forces.  The  experiments  performed  at  the  Naval  Research  Laboratory  and 
Columbia  University  correspond  to  the  regime  covered  by  Case  (3a).  The  e processes  have 
been  listed  in  order  of  increasing  interaction  strength.  Section  IV  contains  a discussion  and 
illustrates  an  example  which  applies  to  generation  of  coherent  visible  radiation.  Section  V gives 
a description  of  the  pertinent  experimental  work,  and  conclusions  regarding  future  work,  both 
experimental  and  theoretical,  are  presented  in  Section  VI. 


S r 


? 


II.  PHYSICAL  MODEL  AND  ANALYSIS 

We  consider  a completely  non-neutralized  electron  beam,  the  motion  of  which  is  treated 
fully  relativistically.  Only  spatial  variations  along  the  z axis  will  be  considered  for  the  electron 
beam,  pump  field  and  scattered  fields.  The  scattered  fields  consist  of  electromagnetic  as  well  as 
electrostatic  waves.  The  amplitude  of  the  static  magnetic  pump  is  taken  to  increase  adiabatically 
in  space,  reaching  a constant  value  for  z > 0,  as  shown  in  figure  1.  The  reason  for  including 
this  feature  in  our  model  will  become  clear  when  the  form  for  the  equilibrium  electron  distribu- 
tion function  is  chosen.  We  choose  as  the  pump  wave  a magnet;c  field  of  the  form 

B„  = B„(ex  cos  k0z  + e,  sin  k„z)  (1) 

where  kn  = 27r//and  / is  the  period  of  the  static  field.  The  pump  amplitude  is  considered  to  be 
a function  of  z for  z < 0 and  constant  for  z > 0.  Although  the  field  in  (1)  does  not  satisfy 
V x B„  = 0,  we  may  consider  it  to  be  a good  approximation  near  the  axis  r = 0 of  the  exact 
helically  symmetric  field  that  is  curl-free  and  varies  as  ljk„r).  If  the  beam  radius  is  taken  to 
be  much  smaller  then  l/2n,  the  actual  field  can  be  well  represented  by  (1).  In  our  model  the 
self  electric  and  magnetic  fields  of  the  beam  are  neglected. 

a)  General  Formalism 

The  pump  field  in  Eq.  (1)  is  derived  from  the  vector  potential 

A„  = - B„M„  (2) 

where  B„  = V x A„.  The  steady  state  equations  of  motion  for  the  relativistic  beam  particles  in 
the  field  of  Eq.  (2)  are 

dpx  fbz  - bAjbz, 

c 

dpjdz  — dAm./dz.  (3) 
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p:  dp/d: 


_ M 


Pv 


34, „ 

82 


+ p, 


8/j„, 

02 


Equations  (3)  yield  the  following  three  constants  of  the  motion  for  the  individual  beam  parti- 
cles: 


<*(.px,z)  = Px  ~ Aox(z)  - px  + - | cos  k0z\ 

c ck0 

P(pv,z)  - px  - Aov(z)  - p,  + ■ | sin  *„z;  (4) 

C 

u(p)  - (px2  + p2  + p?) 1 /2; 

where  a and  /3  are  the  canonical  momenta  in  the  x and  y directions  and  w is  the  magnitude  of 
the  total  momentum.  These  constants  of  the  motion  will  be  used  later  to  construct  the  equili- 
brium distribution  function.  The  electric  and  magnetic  fields  of  the  scattered  waves  are  denoted 
by  E(z.r)  and  B(z,r)  and  are  given  in  terms  of  the  potentials  4>(z,t)  and  A (z,r)  = {Ax.  Ax.  0) 
by 


E (z,r) 


d<f>  . 1_  8A 

8z  c dt  ’ 


B(z,r)  = V x A. 


(5) 


We  shall  perform  a perturbation  Vlasov  analysis  by  expanding  the  electron  distribution  function 
to  first  order  in  the  potentials  <t>  and  A,  about  the  exact  equilibrium.  The  relativistic  form  for 
the  Vlasov  equation  is 

+ v.  - M(E  + v x (B„  + B)/c)  ■ = 0. 

ot  oz  dp 

The  electron  distribution  function  is  written  as 

/(p,z,r)  =/<0)(p,z)  + /(l)(p,z,r),  (6) 

where  /<0)  and  /(l)  are  solutions  of  the  equilibrium  and  first-order-perturbed  Vlasov  equations, 
respectively: 

v.  8/(0,/3z  - ^vxB/  8/(0)/8p  - 0;  (7a) 

c 

0/(l,/8f  + v.a/(l,/3z  - M V X B0  • 8/'<l,/3p 

C 

- k|(E  + v X B/c)  • 8/(O70p;  (7b) 
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where  p = ymay  and  y - (1  4-  p2/m2c2)'n. 


It  proves  convenient  at  this  point  to  make  a transformation  from  the  independent  variables 
(px,py,p.,z,t)  to  new  independent  variables  (a,  p,  u,  z,  /).  Thus  we  write  the  equilibrium  and 
first  order  components  of  the  electron  distribution  function  as 

gm{a,p,u)  - /(0>(p,2). 


g(1)(«,  p,  u,  z,  t)  - 2.  /)• 


In  terms  of  these  new  independent  variables,  Eq  (7b)  takes  the  form 


bgw/bt  + — bgw/bz  - F L g(0) 
ym0 


where 


Y(a,p,u,z,t)  - \e  I E + 


L ia,p,u,z)-ixJ-  + i]'±  + lJ- 


and  the  dependent  variables  are 


px{a,  z)  = a + Aox  ( 2 ), 

c 

py(p,  z)  = P + Aoy  (z), 

p.(a,p,  u,  z ) = [u2  - p2  (a,z)  - p2  (p.z)] 1/2 
y(tt)  = (1  + u2lm2c2Vn. 


. 


Note  that  the  positive  branch  of  the  dependent  variable  p2{a,p,u,z ) will  be  chosen  in  order  to 
represent  a beam  propagating  in  the  positive  z direction.  The  advantage  of  the  transformation 
in  Eq.  (4)  is  clearly  seen  in  Eq.  (9),  in  which  the  Lorentz  force  term,  correspo  iding  to 
y x B0  • d/’Vdp  in  Eq.  (7b),  has  been  transformed  away.  We  also  exploit  the  symmetry  of  the 
model  by  transforming  to  the  vector  basis  (e+,  e_,  e2),  which  is  given  in  terms  of  (ex,  ev,  e.) 
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(10) 


*±  “ y (*.v  ± ' *,■)• 


Noting  Eq.  (5)  we  may  now  write  the  scattered  fields  E and  B in  the  interaction  region  as 

E(z.  i)  - 

a,  \ I dA+(z)  A d.4_(z)  I 

B(z,  r)  =-  - / — e+ eAe~"“'  + c.c.. 


1 d<£(z)  a,  , . . . 

y e-  + A+(z)e+  + /—  AJz)e_ 


e-'-'+c.c.. 


dz  ‘+  dz  '-r  (11) 

where  *(z.t)  - j*(z)*-'“'+  c.c.,  A±(z,r)  - A ±(z)e±e-'“' + c.c.  and  c.c.  denotes  the  com- 
plex conjugate.  The  chosen  time  dependence  of  the  fields  in  Eqs.  (11)  is  appropriate  after  all 
the  temporal  transients  have  decayed  away  and  the  spatial  dependence  along  the  z axis  has  been 
left  arbitrary.  In  this  new  basis  representation  the  force  vector  F and  the  vector  operator  L 
introduced  in  Eq.  (9)  take  the  form 


F(a, (},u,z,t)  - F+(a, p,u,z)e"“'  e+  + FAa,p,u,z)e"“'e_ 


+ Fz(.a,f3,u,z)e"“'e:  + c.c. 


(12a) 


and 


L (a,0,u,z)  L+(a,p,u,z)e+  + Z._(a,/8,«,r)e_  + Lz(a,  0,u,z)e., 


(12b) 


where 


= M 


i(o  — v:(a,  f3,u,z)  — 
dz 


A ±(z). 


F _ M 
2 


d<Mz) 


1 


dz  y(u)m0c 


, _ . dA+(z)  dAjz ) 

p_(a,0,z)  — — + p+(a.p,z) 


dz 


dz 


L±  = 


d -j-  ,_d_  P±(a,/3,z)  jj_ 
da  df3  u du 


P:(a,  P,u,z)  g 


u du' 


vz(a,fi,u,z ) = pz  (a,  fi,u,z)/ (y (w) m0) , 


and 
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i 


(13) 


P±  ia.) 3,z)  - px  T ip  - a?  /£ -7  - e+'*°\ 

ck„ 


The  perturbed  part  of  the  distribution  function  may  also  be  written  in  the  form 
£<n(a,  0,  u,  z,  t)  - g(i)  (a,  0,  u,  z)e~,wl  + c.c. 


The  Vlasov  equation  for  £(1>  now  becomes 


(- ito  + v.(a,p,u,z)d/dz)  g(l)  e + c.c. 


F • L g<0)  — Hia,fi,u,z)e  4-  c.c.. 


where 


L 2 _L_,  J-L(O) 


+ \\i,0-V:t\A-\  I*'  Wj* 

+ |/<o(p_,4  + + p+A J - cp.-^-l  — — , 

( 3z  J u bu 

The  general  solution  of  Eq.  (14)  is 

f(l>  ia,  /3,  u,z)  = dz'  Mia,  f3,u,z,z')  H (a,  0,  u,  z'),  (15) 

where  we  have  assumed  that  the  beam  entered  the  interaction  region  unperturbed,  i.e., 
g(l>  (a,  p,  u,  0)  - 0,  and  where 


Mia,  /J,  u,  z,  z')  — — - — — exp(/&)r(o,  fi,u,z,z')) 

v.ia.p  ,u,z') 


ria.p.u.z.z’)  = f — ■— izo — . 

J - v.ia.p.u.z) 

Rearranging  Eq.  (15),  we  may  write  the  Fourier  coefficient  of  the  perturbed  distribution  func- 


tion as 
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P u'z)  “ G+(«.P.u,z)  |-£-  + i -^)  + GAa.p.u.z) 

+ G.(a.p.u.z)  -j-  g™(a,p,u)  (16) 

au 

where 

(j±  - f dz'  M(a,p,u,z,z‘)  |/a>  - v.(a,p,u,z')  -£-?  A±(z’),  (17a) 

2 c •'O  I or 

~ f <fe’  M(a,P,u,z,z‘)  f i<o pAa.p.z')  A+(z') 

2c  u l 

+ iu>p+(a,p,z')  A-(z')  ~ cp.(a,p,u,z')  ^ ^ . (17b) 

or 

Careful  inspection  of  the  definition  of  G±  shows  that  the  integrand  in  (17a)  is  a perfect 
differential;  thus,  we  can  perform  the  integration  over  z'  and  obtain 


G ±(a,p,u,z)  = p ± (z)  - A ±(0)  exp(/o>  t (a,/3,w,z.O))  J . (18) 

The  expression  in  Eq.  (15)  for  determines  the  first  order  perturbation  of  the  electron  distri- 
bution by  the  scattered  fields  with  arbitrary  axial  space  dependence,  correct  to  all  orders  in  the 
pump  amplitude. 


The  perturbing  driving  current  density  J(z,/)  wfiich  excites  the  scattered  fields  can  be 
found  from  #(I>  and  is  written  in  general  as 


where 


J(z,r)  = lJ+(z)  e+  + J_(z)e_  -I-  J.(z)  c:]  e + c.c.. 


Ld  f“  -7—  I Pf('a'f’Z\\  g(l\a,p,u,z ) r dadpdu  (20) 

l(z)  J m0  Jo  y (u)  | p.(a,p,u,z)  j p.(a,p,u,z) 


and  we  have  used  the  fact  that 


dlp  — — - — 7 - dadpdu. 

p.(a,  p.u.z) 

The  general  framework  of  the  formalism  is  completed  by  relating  self-consistently  the  scattered 
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fields  in  Eqs.  (11)  to  the  driving  current  density  given  by  Eq.  (19).  This  is  done  through  the 
wave  equations  for  <f>(z,t)  and  A±(z,/)  which  are  given  by 


■ 


[ 


! 


d2<t> 

dtdz 

J _ jP_ 

c2  dr2 


4tt  J., 


J± 


(21a, b) 


It  should  be  mentioned  at  this  point  that  the  terms  G±  which  contribute  to  gU)  and  thus 
to  J do  not  appear  in  any  previous  analyses  on  FELs  to  our  knowledge.  This  is  due  to  the  fact 
the  previous  analyses  do  not  use  the  exact  equilibrium  distribution  function  in  their  treatment 
of  the  problem. 


As  already  mentioned,  the  ponderomotive  potential  wave  which  results  from  the  coupling 
of  the  pump  and  scattered  wave  is  responsible  for  the  longitudinal  bunching  of  the  beam  elec- 
trons. The  modulated  beam  together  with  the  pump  field  produces  a current  density  which 
drives  the  scattered  wave.  The  most  straightforward  way  to  appreciate  the  bunching  process 
associated  with  the  ponderomotive  potential  is  to  consider  the  form  of  the  Hamiltonian  for  a 
particle  in  the  combined  fields  of  the  pump  and  scattered  waves.  The  Hamiltonian  of  such  a 
particle  with  zero  transverse  canonical  momentum  is  given  by 

H = /c2p?  + |e|2  (A„(z)  + \(z,t))2  + m2c*  - \e\<t>  ( z,t ) 

where  A„(z),  A(z,r)  and  $(z./)  are  found  from  Eqs.  ( 2 ) and  (5).  Since  |A(z,r)|  <<  |A„(z)| 
we  can  expand  the  Hamiltonian  to  first  order  in  A(z,r)  and  write  it  in  the  form 


H = y/c2p?  + \e\2kHz)  + m2c*  - |e|(<Mz,r)  +<*w(r.')J 
where  <t> ^(z.t)  is  defined  as  the  ponderomotive  potential,  given  by 

|e|A„(z)  • A(z,r) 


•y^oC 


and 
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■pn 
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y = y/\  +p:2/m,}c2+  \e\2A*/(m,}c*). 


In  this  expanded  form  we  see  that  the  transverse  component  of  the  scattered  field  given  by 
A (z.t),  when  coupled  to  the  pump  field,  plays  a role  similar  to  that  of  the  space  charge  potential 
<t>(z,t).  The  ponderomotive  potential  <£,„„„/(-.'),  has  the  effect  of  axially  bunching  the  beam. 
This  bunching  process  in  turn  gives  rise  to  the  space  charge  scalar  potential  4>(z,i).  Hence,  the 
ponderomotive  potential  drives  the  space  charge  potential.  For  future  reference,  we  note  that 
the  ponderomotive  potential  can  be  put  into  the  form 

</w(z-'>  = ^ + cc-  (22a) 


where 


and  n„=  \e\BJm„c. 


,,  . . 
<t>Pomi(z)  = e " A+(z)  + e A Az) 


(b)  Evaluation  of  the  Current  Density 


In  order  to  obtain  an  explicit  expression  for  the  driving  current  density  given  in  Eq.  (19) 
we  need  to  specify  a form  for  the  equilibrium  distribution  function.  As  was  mentioned  earlier, 
the  pump  magnetic  field  amplitude  B„(z ) builds  up  adiabatically  and  reaches  a constant  value 
inside  the  interaction  region.  We  assume  that  the  electron  beam  is  generated  far  to  the  left  of 
the  interaction  region,  where  the  pump  field  vanishes.  In  general  the  beam  will  be  produced 
with  a thermal  momentum  spread  in  the  transverse  as  well  as  the  longitudinal  direction.  How- 
ever, the  transverse  momentum  spread,  PLlll>  may  be  neglected  if  it  is  much  less  than  the 
ordered  transverse  momentum  \e\B„/ck,„  induced  by  the  pump  field.  Unless  the  pump  ampli- 
tude is  unusually  small,  the  condition,  Puh  « \e\B,Jck„  can  be  easily  satisfied.  Then  the 
transverse  momentum  is  proportional  only  to  the  pump  amplitude  and  the  solutions  of  the 
equations  of  motion  are  accurately  given  by  a - 0 = 0 in  Eqs.  (4).  Thus,  for 
P1,h  « |c|  B„/ck,„  we  may  choose  the  equilibrium  distribution  function  to  be  of  the  form 


£<0,(a,  P.  u ) - n„f>(a)  8(0)  g„(u), 
15 


(23) 


where  n„  is  the  ambient  beam  density  far  to  the  left  of  the  interaction  region.  The  normaliza- 
tion of  g„(n)  is  such  that 


r 

•'ll 


du 


u.  (u) 


g„(u)  = 1, 


where  u.(u)  = pr(a  = 0,  p - 0,u,  z)  = (u1  - \e\2B,}/c2k,})'n.  We  may  remark  that  for  a - 
P “ 0,  each  particle  has  constant  axial  momentum.  This  is  a rather  important  property  of  the 
pump  field  ip  Eq.  (1).  Since  the  equilibrium  axial  momentum  and  velocity  do  not  oscillate  as  a 
function  of  z the  pump  field  does  not  introduce  an  effective  longitudinal  spread  in  velocity 
which  would  limit  the  effectiveness  of  the  scattering  process.  The  appropriate  functional  form 
of  £„(w)  will  depend  on  the  scattering  process  being  considered;  we  shall  leave  it  general  for 
the  moment.  Using  the  distribution  in  Eq.  (23)  it  is  straightforward  to  evaluate  the  Fourier 
amplitudes  of  the  driving  current  density  in  Eq.  (20).  The  details  of  this  calculation  are  given 
in  Appendix  1.  Defining 


'fMz)  = 


~ ± k. 


we  find 


JAz) 


47 T 


r 

Jo 


yU.C 


A ±(z)  - A AO)  e 


i<o:/ v. 


, P±  to 

+ T"  T — e 

2 k„v. 

|2 


A +(0)  (l  -e  k':\  - /i_(0)  (l  - e j] 


+ ^ [AAz)e-'**i:' - -^(O)^]) 

r-  ja^(z’)  , uP±  f-  , „ 

Jo  IF"' 


f-  \U(Z')  ^ "h  | . a 


+ A Az')e 


dz'  f- 
6w 


gn  du 

(24a, b) 


f \~  S ')  - '*-<»»  («'v-  >11 


- ™-  X'l4^2 + **  *•£ 


g„(u ) du,( 24c) 
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where  toh  = (4n \e\1njm,) 1/2  is  the  beam  plasma  frequency,  = njyk„v.  is  the  transverse 


beam  velocity  normalized  to  v:,  fl„  - |e|  B„/m„c,  and  v.  - u./ym„. 


It  is  informative  to  consider  at  this  point  the  current  density  J(z,  t ) in  the  absence  of  the 
pump  field  in  order  to  find  the  pump-free  eigenmodes  of  the  system.  Setting  B„  »0  we  find 
from  Eqs.  (24)  that 


- . tt)>,  r°°]  ..  ■ 

J±(z)  I = —— J — A ±(z)  - A ±( 0)  f 1 g„  du, 

a„-o  4i tc  y 1 " 


Uz)  | - + ~ m„  f f 

5-0  8tt  •'«  Jo 


dMz')  d . 

e dz  t;  gn  du- 


where  v.  | = u/ym„.  For  the  moment  we  will  consider  the  eigenmodes  for  z » 0.  In  this 

- o 

case  we  can  take  A ±(z)  and  <Mz)  to  be  plane  waves  of  the  form  /i±(z)  = /i±(0)  exp(/  k±z) 
and  <j>(z)  = 4>(0)  exp(/  k z ) where  k ± and  A are  real  constants.  Neglecting  the  free  streaming 

terms  in  Eqs.  (24),  i.e.,  exp  t zj,  because  they  eventually  phase  mix  away,  we  obtain  for 


z » 0 


/•“  g„  du 


1 1 \ 

J-(z)  = — — m, 

oTT 


9g„/du  du 

Ooj  J — <Mz). 

(o  — v.  k 


Combining  the  current  densities  given  by  Eqs.  (19)  and  (26)  with  the  wave  equations  in  (21), 
we  arrive  at  the  conventional  dispersion  relations  for  the  transverse  electromagnetic  field  and 


the  electrostatic  field: 


D (to,  A ±)  A ±(z)  = 0; 


(1  + x(“,  A))  <j>(z)  - 0; 


where  x is  the  electron  susceptibility 


f ,\  C m«  dZn/Bu  du 

X(<o.  A)  - — I — 

A J to  — v.  k 


i 

! 


and 

D(a».  k±)  = a>2  - c2k2±  - f gjy  du.  (28b) 

It  is  now  possible  to  show  that  for  the  doppler-upshifted  backscattered  wave  we  can  con- 
sider the  A + or  A.  wave  separately.  In  the  presence  of  the  pump  field  it  will  be  shown  that  the 
products  D(w,k±)  /4±(z)  are  no  longer  zero,  but  are  equal  to  a small  coupling  term  propor- 
tional to  B2.  In  order  for  -4+(z)  to  be  excited  it  is  necessary  for  DU,  k+)  = 0;  likewise,  for 
A~(z)  to  be  excited  DU,A;_)  must  be  very  small.  It  can  be  shown  by  carefully  scrutinizing 
Eqs.  (24)  that  k+  and  are  connected  by  the  relationship  k+~  k_  - 2 k„.  Because  of  this 
relationship,  the  functions  D( w.  k+)  and  DU,  kj  cannot  be  simultaneously  zero  for  the  same 
backscattered  upshifted  frequency.  Therefore,  A+  and  A_  are  independently  excited.  Further- 
more, by  simply  changing  the  sign  of  oj  and  either  k+  or  k_  we  see  that  the  behavior  of  the  two 
waves  are  identical.  Hence,  in  the  remainder  of  this  paper  we  shall  arbitrarily  take  the  wave  at 
(o>,  k+)  to  be  resonant,  i.e.,  D(w,  k+)  ~ 0,  and  the  (to,  k_)  wave  to  be  nonresonant,  i.e., 
DU,  k J ^ 0.  This  amounts  to  keeping  only  the  A+  wave  and  neglecting  the  A_  wave.  We 
will  now  evaluate  the  growth  rates  and  saturation  level  for  various  scattering  regimes. 


III.  GAIN  AND  SATURATION  LEVEL  IN  VARIOUS  FEL  REGIMES 


Case  1.  Low  Gain,  Tenuous-Cold-Beam  Limit 


As  our  first  example  we  will  consider  the  low  gain,  (,.<>..  "short  cavity")  scattering  process 
in  which  collective  effects  are  unimportant.  This  limit  corresponds  to  the  parameter  regime  of 
the  recent  experiments  carried  out  at  Stanford  University  with  highly  relativistic  48  MeV), 
low-current  (<  2.4A)  electron  beams.  Collective  effects  are  manifested  by  the  space  charge 
potential  <t> , which  arises  from  the  charge  bunching  under  the  action  of  the  pondermotive 
potential  4^.  In  the  very  tenuous  beam  limit  the  scalar  potential  is  much  less  than  the  pon- 
aermotive  potential.  It  can  be  shown  that  collective  effects  can  be  neglected  when  the  magni- 
tude of  the  electron  susceptibility  x is  small  compared  to  unity.  In  this  limit  we  can  set  <t>  — 0 
and  represent  the  scattered  electromagnetic  wave  in  the  form  (see  Eq.  (ID) 

A + (z)  = /i+(0)  e'  *+<" 1 d'  (29) 

where  k+(z)  = k0+  + 5&(z),  k0+  » |SA | is  real  and  constant  and  Bk  is  complex  and  a slowly 
varying  function  of  z.  Using  this  form  for  the  field  we  find  from  Eq.  (24a), 


•7+(z>  = — 


/ —I  1 + (l  - _ BL  _2»_  (j  _ -'V) 

yu:  I 2 l ) 2 k0  v:  tie; 


• » ^ z ) + k z ) /*'  , i u z ) + k z ) f)  » . C l si-  • 

+ >~2  Jn  dz  e • -2-  go  m+(0)  e du.  (30) 

where  i \>(u,  z ) - ^ k+(z')  dz'  - cjz/v..  In  the  context  of  Eqs.  (29)  and  (30)  the  wave  equa- 
tion (21b)  can  be  used  to  derive  the  following  dispersion  relation: 


k0\  - ~ + 2k0+Bk(z)  = n(u>,  k0+.  z). 


where 


d 


1 „)  (i  — e 

— PS  mnv-m  


e ) a . 

Z T~  du 

oj  — v.k  du 


and  where  k = k„+  + k„. 


Note  that  in  obtaining  (31)  we  used  the  facts  that  in  the  present  scattering  limit, 
| SA:  (z)  | <<  k0+,  Jq  8k(z')dz'  <<  1 and  81c  (z)  is  a very  "weak"  function  of  z,  so  that  z) 
has  been  approximated  in  the  RHS  of  Eq.  (32)  by  (/c0+  — <o/v.)z.  Solving  for  the  imaginary 
part  of  Eq.  (31),  we  find 


Im(Slc)  = - ^ Jf  y P±m0o 


sin  ( ou/v . — k)z  3 
(o >/v.  — k)  du 


~ yy  |(1  + /3i2/2)sin(a»/v-  - k„+)z 
+ y-  [sinGu/v-  - k„+)z  - sin(a>/v.  - lc)z||L„</w. 


The  total  gain  in  amplitude  of  the  scattered  wave,  over  the  interaction  length  L,  is  given  by 
G,  = - JqL  lm  (8k  (z'))dz‘ 

= r°°  El  2 sin(/c  - wo  l/2  2 _a_ 

2 k0+  Jo  4 mo<1>  (k  - a>/ v.) L/2  ~du 


u Pi  oj  L sin2(lc  -(o/v-)L/2 

y u.  2 v.kn  (k  — (o/  v.)  L /2 


. , Pi  W sin2(k0+  — ai/v.)L/2 

2 v.k„  (k0+  — <xi/v.)L/2  gndu‘ 


Equation  (34)  is  a complicated  expression  for  GL  which  depends  on  gju).  If,  however,  g„  is 
such  that  the  phase  of  the  sine  functions  in  Eq.  (34)  does  not  vary  appreciably  over  the  range 
of  v-  for  which  g„  is  non-zero,  then  to  evaluate  the  integral  we  may  take 


g„(u)  = — 8(u  - u„). 
u 


The  condition  for  the  applicability  of  (35)  is  that  the  axial  thermal  velocity  spread  V„,  of  the 
beam  satisfy 
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v>h  <<  ~TTJ  “ ~T~-  06) 

y Lk0L  L 

The  thermal  energy  spread  E„,  of  the  beam  in  general  is  given  by 

£)/,  = y,,y  ?om„V;0Vlh,  (37) 

and  for  a highly  relativistic  beam,  inequality  (36)  may  be  rewrite  as 

~ « >1  (38) 

L 

where  E„  — (y„  - l)m„c2  = y„m„c2  is  the  kinetic  energy  of  the  particles.  In  the  experiments 
to  date  this  condition  is  satisfied  by  roughly  one  order  of  magnitude. 


For  a highly  relativistic  beam,  we  may  approximate  v.  = c,  u.  = w,  klt+  = 2 y2k„.  Then 


using  Eq.  (35),  the  total  gain  becomes. 


Gl  - t4-  k„L 
*-y:o 


,,  r,2  3 sin6>0  2 , sin20o 

2 3 90  0O  +2 


2 2,  sin2(e„  + k0L/2) 


- 0 + Po±y& 


(0O  + k0L/2) 


2y.„  ) k°L  Fi9°’  k°L'  Po±  ^ 

where  f - (ob/(J~yock0)  is  the  beam  strength  parameter,  y0  • y(u0),  y-0  = y.(u„),  0^ 
= ft.(w0)>  and  0„  = (<*>/ v;rt  — Ar„+  — k0)L/2.  In  deriving  Eq.  (39)  we  have  used  the  fact  that 

3 _ _ (o  L u 3_ 

v?  2 ' m0y  u.  y2  30, , 

In  Figures  (2a,b)  we  exhibit  the  function  F(90,  k0L,  p01  y:o)  for  the  cases  where  L = 160/ 
and  L = /.  If  the  interaction  region  consists  of  a large  number  of  pump  periods,  k„L  >>  1, 
the  first  term  of  the  gain  expression  in  Eq.  (39)  is  dominant.  For  this  case 


* *4*  «■«*£ 


The  function  3/3 9n  (sin  0o/0„)2  has  a maximum  value  of  0.54  when  0„  = -1.3,  therefore,  the 
maximum  value  of  gain  is 


G,  ~ WVL  «oL)\ 
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These  expressions  for  the  gain  are  similar  to  results  found  previously  by  Motz  (1951), 
Sukhatme  and  Wolff  (1973),  Colson  (1976),  Hopf  et  al.  (1976a),  and  Kroll  and  McMullin 
(1978). 


A detailed  calculation  of  the  efficiency  in  the  free-eleclion-laser  low-gain  regime  was 
presented  by  Hopf  et  al,  who  found  that  the  efficiency  was  approximately  given  by 

Tj  = (Xj/L)  y2  (42) 

where  Xv  -*  lit lk+.  It  is  possible  to  obtain  this  result  from  general  arguments  based  on  elec- 
tron trapping  in  the  longitudinal  wave  associated  with  the  electron  beam.  (The  total  longitudi- 
nal wave  in  general  consists  of  both  the  space  charge  potential  and  ponderomotive  potential 
waves.)  The  following  discussion,  which  will  be  used  to  estimate  saturation  efficiencies  as  well 
as  the  saturated  electromagnetic  field  amplitude,  may  be  applied  to  all  the  cases  involving 
scattering  from  a cold  beam  in  which  trapping  saturates  the  process. 

For  maximum  spatial  gain,  the  phase  velocity  vrh  = a >/Re(k)  of  the  longitudinal  wave  and 
the  directed  axial  beam  velocity  v.„  are  related,  in  the  linear  development  of  the  instability,  by 
an  expression  of  the  form 

V,  - v.„  = -Av  (43) 

where  Av  > 0 is  a small  term  which  depends  on  the  particular  scattering  regime  being  con- 
sidered. In  the  low-gain,  tenuous-beam  limit  we  have  lust  discussed,  for  example,  the  gain  GL 
maximizes  at  0„  = - 1.3;  this  condition  may  be  writ  en  in  the  form  (43).  In  the  high-gain 
regimes  to  be  discussed  later,  in  which  the  scattering  process  may  be  described  by  a dispersion 
relation  homogeneous  in  space,  a condition  of  the  form  (43)  will  be  seen  to  describe  the  condi- 
tion of  maximum  coupling  between  the  beam  and  the  backscattered  wave.  Since  wave  growth 
will  occur  for  Av  > 0,  we  can  interpret  this  condition  as  stating  that  the  phase  velocity  of  the 
longitudinal  wave  must  be  less  than  the  axial  beam  velocity  in  order  for  beam  kinetic  energy  to 

be  converted  into  wave  energy.  The  wave  amplitude  increases,  at  the  expense  of  the  beam 
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kinetic  energy,  until  the  beam  electrons  become  trapped  at  the  bottoms  of  the  potential  wells  of 
the  longitudinal  wave.  At  this  point  the  wave  reaches  its  maximum  amplitude  and  the  average 
beam  kinetic  energy  reaches  a minimum.  The  average  axial  beam  speed  at  saturation,  v.s,  is 
approximately  given  by  the  trapping  condition 

V - v-v  = + Av  (44) 

Therefore,  at  saturation  the  kinetic  energy  lost  by  the  beam  corresponds  to  a decrease  in  axial 
beam  speed  of  v;„  - v>.  s = 2 A v.  It  should  be  noted  that  when  the  beam  electrons  are  trapped  at 
the  bottom  of  the  potential  wells  of  the  wave  they  are  modulated  spatially  but  are  fairly 
monoenergetic.  The  decrease  in  beam  energy  can  be  equated  to  the  increase  in  wave  energy. 
In  this  way  the  efficiency  and  saturated  field  amplitude  of  the  scattering  process  can  be 
estimated. 


In  the  case  of  a relativistic  electron  beam  the  change  in  kinetic  energy  at  saturation,  when 
the  beam  particles  are  deeply  trapped,  is 
A^K.E.  = (y(v.„)  ~ y(v.j))  m J.C2 


„ 6yU  = vM)  . , _ 

— 2 Av  mnc  = 2y„y.„  m0v.o Av. 


(45) 


The  efficiency  is  then  given  by 


11  = 7 — ~nT — 7 — ^v/c 
(y„  - l)  mo? 

where  the  last  expression  in  (46)  is  for  a highly  relativistic  beam,  v;o  ~ c.  The  amplitude  of 
the  scattered  electromagnetic  fields  at  saturation  can  now  be  obtained  by  the  following  argu- 
ment. 

In  the  steady  state,  conservation  of  total  energy  flux  states  that 


JL 

dz 


nv.( y - 1)  m„c2#.  + E x B 

4ir 


0. 


Integrating  from  the  input  of  the  interaction  region,  z “ 0,  to  the  point  where  the  scattered 
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field  is  maximum,  z - zw„  gives 

I'M*  = zM,)|2  = |/i+(z  = 0)|2  + -n~C-  A £ke 
~ |/i+(2  “0)!2+ 


(47) 


In  obtaining  \A+(z  = zM,)|  we  used  the  fact  that  the  fields  were  given  by  Eqs.  (11),  A_  - 0, 

E x B ~ (<o2/c2)  M+(z)|2  and  v.0  - v.,  = c.  The  magnitude  of  the  transverse  electric  field  at 
saturation  becomes 


|E| 


I'M*  = zM,)|. 


(48) 


The  results  for  the  efficiency  and  saturated  transverse  electric  field  in  Eqs.  (46)  and  (48) 
apply  to  those  scattering  processes  in  which  electron  trapping  is  responsible  for  saturation.  To 
apply  these  results  to  the  low  gain  scattering  process  in  Case  1,  we  note  that  the  condition  for 
maximum  gain  over  a distance  zM,  - L is  given  by  90  - - 1.3,  which  when  put  in  the  form  of 
Eq.  (43)  becomes 


Vph  ~ V--°  “ ~kL 

Substituting  Eq.  (49)  into  Eq.  (46)  the  estimated  value  of  efficiency  is 

t)  ~ 5.2  y?J(kL)  ~ y?0  ks/L, 


-2.6 


v.„ 


- Av. 


(49) 


(50) 


in  agreement  with  Eq.  (42),  where  X,  = 2n/k0+  ~ 2 n/k  is  the  wavelength  of  the  scattered 
wave.  The  assumption  implicit  in  the  above  argument  is  that  the  length  L and  the  energy  lost 
by  the  electrons  are  matched  so  that  the  trapping  occurs  at  the  end  of  the  system.  The  energy 
lost  by  the  electrons,  however,  goes  into  the  scattered  wave  and  depends  upon  the  amplitude  of 
the  input  signal.  To  find  the  field  amplitude  necessary  to  effectuate  this  energy  change  in  the 
length  L we  note  that  in  the  low  gain  limit  GL  « 1,  we  have  \A+(z  - zM,  - T)| 
~ | A+(z  = 0)|  (1  + Gl)  and  thus  Eq.  (48)  gives 


M\(r  -0)|2 


2nn„ 


U/c)2  (G, 

At  maximum  gain  conditions  we  have,  using  Eq.  (41), 


-2L-I  y„m„c2. 
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M+(z  -0)|  - >/5T 


Yo 

m„c 2 

1 1 

Y:a 

kl 

Poi  ( k„L )2 

(51) 


Since  M+(z  -zM,)|  — \A+(z=  o)|,  the  saturated  transverse  electric  field  amplitude  can  be 
estimated  from  Eqs.  (48)  and  (51)  as 


|EL,  ~2yfS2 


m„c‘ 


kl 


Yo 


" Poi(k„L)2 


(52) 


Equation  (52)  gives  the  required  input  signal  necessary  to  achieve  the  efficiency  in  Eq.  (46)  at 
maximum  gain  conditions.  This  input  signal  level  is  required  for  optimum  efficiency.  If  the 
input  signal  is  weaker  than  this  value,  saturation  of  the  wave  will  not  occur  within  the  length  L. 
If  it  is  stronger,  saturation  will  occur  before  the  beam  reaches  the  end  of  the  interaction  region, 
resulting  in  smaller  total  gain. 


Case  2.  High  Gain,  Tenuous-Thermal-Beam  Limit 


We  now  consider  the  high  gain  limit  with  a tenuous  thermal  beam;  this  limit  is  the  rela- 
tivistic analog  of  stimulated  Compton  scattering.  We  again  neglect  the  space  charge  potential 
but  now  assume  that  the  wave  e-folds  many  times,  hence  |/4(z)|  » |/4(0)|  except  near  z — 
0.  Since  for  this  case  -Im  \ k+(z')  dz'  » 1,  the  driving  current  is  found  from  (30)  to  take 
the  form 


J+(z)  - — 


«•»  r°° 

u 

i + El 

Pi  d/du 

4irc  •'» 

y“= 

1+  2 

(i)  v:.m0 

2 c o — v.  k 

ik 


g„  du  /4+(0)  e *'  (53) 

where  k+  is  complex  and  independent  of  z.  Substituting  Eq.  (53)  into  Eq.  (21b)  gives  the 
dispersion  relation 


where  k = k+  + k„  = kr  + ik,.  In  the  Compton  scattering  regime  the  phase  velocity  of  the 
ponderomotive  wave  is  resonant  with  a small  fraction  of  the  thermal  beam  electrons  (see  Fig. 


3).  We  now  consider  the  case  where  the  thermal  width  of  the  particle  distribution  function  is 
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large  compared  to  the  width  of  the  wave  distribution  (w-v./r)-1.  If  the  thermal  velocity  spread 
of  the  relativistic  beam  is  such  that 

Klh  » | c kjk, |.  (55) 

then  the  pump  wave  scatters  from  single  particles.  When  the  inequality  (SS)  is  satisfied  we  can 

take  the  imaginary  parts  of  both  sides  of  Eq.  (54)  and  use  the  fact  that 


Im 


1 


c o — v.k 


l*rl 


8(v.  - v/kr) 


as  k,  approaches  zero  from  negative  values.  Because  v.  = c(w2  - Uq1)u1/(u2  + m2c2)U2,  we 
can  write 


x yy;m„u: 

8(v.  - to/k, ) = 8(w  - uph ) 


(56) 


where  dvjdu  « u/(y  y?m0u.),  u > 0,  u,„  = yph{Pl,  m,l  c2  + u0i)'/2,  u.  = (i/2  - u2±)U2y 
Pph  = <o/k,c,  ypt,  = (1  - P!h)-'/2.  and  u^  = (|e|  Bjck0 )2.  The  imaginary  part  of  k,  from  Eq. 
(54)  using  Eq.  (56),  is  given  by 


7T 


P±  V;  m„  0) 


y y?  m0  u.  dg„ 


(57) 


4 k,c2  |kf|  '"'i  '■  u du 

For  a drifting  electron  beam  with  average  total  momentum  we  can  take  g0(u)  to  be 

Maxwellian-like  of  the  form 


„ (W)  = UAU  ..ttUI/J) 

8o(u)  ^ u,h  e 


(58) 


such  that  s:  du  g0  u/u.  = 1.  Substituting  (58)  into  (57)  we  get  for  a highly  relativistic  beam 
the  result 


r a2  y mo  _v2 

**r  ~ur 


(59) 


V 


where  x = («  — w„)/(>/2  Ulh).  The  function  xe~%1  has  a minimum  value  of  -0.43  when 
x - -1/727  i e.,  upll  = u„  - Ulh  Since  for  a highly  relativistic  beam  the  thermal  energy  spread 
£,/,  is  very  nearly  equal  to  cU„,y  the  maximum  spatial  growth  rate  becomes 

. 12 


l*,L.,  = — 


E„, 


(60) 


where,  as  before,  E„  = y„mc 2,  and  f - a >J(yJy'o  ck„). 
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These  results  for  the  growth  rate  agree  with  those  of  Lin  and  Dawson  (1975),  Hasegawa 
et  al.  (1976),  Kroll  and  McMullin  (1978),  and  Sprangle  and  Drobot  (1978). 

In  order  to  estimate  the  efficiency  in  this  stimulated  Compton  regime,  we  first  note  that  only  a 
small  fraction  of  the  beam  particles  are  resonant  with  the  unstable  wave.  This  fraction  is  deter- 
mined by  the  width  of  the  spectrum  of  the  spatially  growing  wave.  The  width  of  the  wave  spec- 
trum in  velocity  space  is  given  to  a good  approximation  by 


8v,  = \2ck,/k\  « ylh. 

In  momentum  space  (see  Fig.  4)  this  width  becomes 


(61) 


6m, 


__  du 
dv. 


8v»  ~ |2y0  yf0  c k,/k\ 

~|y„  m„  c kjk„\  « U„,  (62) 

As  the  backscattered  electromagnetic  wave  grows,  a small  fraction  of  the  beam  particles,  those 
with  momenta  between  uph  - huw/ 2 and  unh  + 8m„,/2,  become  trapped  in  the  potential  associ- 
ated with  the  ponderomotive  wave.  As  this  happens,  the  distribution  function  begins  to  flatten 
out  near  u = uph  and  dg/du  decreases  until  it  vanishes.  A plateau  is  then  formed  and  the  dis- 
tribution becomes  stable  with  respect  to  the  initially  unstable  wave  because  dg/du  | - 0 

« - v 

(see  Fig.  4).  The  change  in  kinetic  energy  between  the  initial  and  final  distribution  functions  is 
given  by 

where  we  have  used  the  fact  that  unh  = u„  — U,/,  for  the  maximum  growth  rate.  Since 
8 uw  « Ulh,  evaluation  of  Eq.  (63)  gives 


8 E, 


K E. 


~ n°  e-\n  \ 8u* 
y/2^  ‘ \2Ulh 
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du 
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Wt 
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where  dujU,h  is  evaluated  from  Eq.  (62)  with  k,  given  by  Eq.  (60)  and  is  given  by 
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(64) 


(65) 


*"JUm  - e P liEjEj'  « 1. 


The  efficiency  using  Eq.  (64)  is  given  by 


1 E,h 


n„(y0  - l)  mBc2  4>/2tt  ( U,h  J E„  ' 

Using  the  expression  for  the  efficiency  in  (66)  and  Eq.  (47)  together  with  (48),  we  find  that  the 
electric  field  amplitude  at  saturation  is 

Irl  m°c2  f i ® Uw  SJEg 

|E'“  - TT r* s *•  a:  <6,) 

The  requirement  on  the  beam  thermal  velocity  is  stated  in  Eq.  (55).  In  terms  of  the  ratio 
EJE„ , this  condition  becomes 


EJE„  » £ 2 pi 

where  Eqs.  (55),  (60)  and  (37)  have  been  used. 


Case  3.  High  Gain,  Dense-Cold- Beam  Limit 

In  this  limit,  we  shall  retain  the  space  charge  scalar  potential  <f>  which  arises  from  charge 
bunching  by  the  ponderomotive  potential.  Because  the  scattered  fields  are  assumed  to  e-fold 
many  times  within  the  interaction  region  we  may  neglect  boundary  terms  at  z - 0 and  choose 
0 and  A + to  have  the  form  (see  Eq.  (11)) 


4>(z)  = <£(0)  , 

A+(z)  - A+(0)e'k+ 


where  Ar+  and  k are  complex  constants. 


Substituting  the  fields  in  (69)  into  Eqs  (24)  and  taking  the  high  gain  limit,  i.e., 
|exp(/  k z)\  » 1 and  |exp(/  k + z)\  » 1,  the  driving  current  densities  become 

J"  f(-^-  (2  + el)  - - f ^,(0)  ,*-• 

8ff  Jo  yu.c  O)  - v.(k+  + k„)  du 

P^nV.kd/du  „ ,(*_*)- 

- -~lv_k  *<°>  e " *.*•.  (70a) 
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du 


du. 


(70b) 


I w — p-Jr  co  — v.(k+  + /c0) 

We  see  immediately  from  Eqs.  (70)  that  the  selection  rule  between  the  wavenumbers  of  the 
three  waves  is 


* - k+  + k» • (71) 

The  coefficients  of  the  driving  currents  3+  and  3:,  using  Eq.  (70),  can  be  written  in  the  form 
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where 
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Substituting  Eqs.  (69)  into  the  wave  equations  shown  in  (21)  gives 

<M0)  eikz  = i+(0)  e‘k:, 

ck„  1 + *(oj  ,k) 


n(  u \ ft° 


S(oi,k) 


'MO)  e'k+:  = - ~ Xy(oi.k)  «(0)  e'k*:  (73) 
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The  first  of  the  equations  in  (73)  can  be  put  into  an  illuminating  form  by  noting  that  to  a good 

approximation  y can  be  taken  out  of  the  integral  in  the  definition  of  *y((u,  *),  since  it  is  a 

slowly  varying  function  of  u.  With  this  approximation  together  with  the  definition  of  the 
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ponderomotive  potential  in  Eq  (22)  we  find  that 


0b)  - -X(w-  k)/2  4 ,(2) 

1 + x( <o.  k) 


(74) 


where  <6(z)  — 0(0)  exp  ( ikz ).  In  Eq.  (74)  we  see  explicitly  the  tact  that  the  ponderomotive 
potential  drives  the  space  charge  potential. 

Combining  the  two  equations  in  (73)  we  arrive  at  the  fully  relativistic  high  gain  dispersion  rela- 
tion 
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where 


Q(<o,k)  - S(a),k)  + 

to2  1 + x(o).  k) 

Again,  we  take  a monoenergetic  electron  beam  represented  by  the  distribution  function 


8n  = 6(m  - u„) 

where  u„  is  the  total  beam  momentum.  For  such  a beam  we  find  that 
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D(<o,k)  — <d2- c2k2-a>2/y„  , 

and  the  dispersion  relation  becomes 
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where  y0  - y(u„),  v.„  - v.(u„),  and  y:o  - y 


Note  that  in  the  absense  of  the  pump,  fl„  — 0,  we  recover  the  uncoupled  dispersion  relations 
(27)  for  both  the  electromagnetic  and  electrostatic  waves.  In  the  presence  of  the  pump,  how- 
ever, these  two  modes  are  coupled.  Since  the  electromagnetic  waves  in  the  presence  of  the 
pump  field  approximately  satisfy  the  pump  free  dispersion  relation  D(<o,  k+)  = 0,  we  can 
replace  £)(&>.  k ) on  the  right  hand  side  of  Eq.  (77)  by  —2kk„c2.  The  dispersion  relation  in 
(77)  may  now  be  re-cast  into  the  form 
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We  may  now  discern  two  limits  for  the  solution  of  (78),  depending  upon  the  pump  strength. 


(a)  W 'eak-pump  regime 


When  the  pump  amplitude  is  smaller  than  a critical  value,  which  we  shall  determine 

below,  we  may  take  the  interacting  longitudinal  and  transverse  waves  to  be  approximately  nor- . 
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mal  modes  of  the  system.  The  positive-energy,  fast  space-charge  mode  of  the  beam  is 
described  by  k = o>/v-„  — k,  while  the  negative-energy  slow  space-charge  wave  has 
k = <t)/ vzo  + k.  It  is  the  interaction  of  the  slow  wave  with  the  pump  that  leads  to  instability  of 
the  backscattered  electromagnetic  waves.  The  coupling  of  the  negative-energy  space-charge 
wave  and  the  transverse  electromagnetic  wave  is  depicted  in  Fig.  5.  Thus,  in  the  case  of  a weak 
pump  we  can  set 


k — kn  -4-  K + bk  — oj/v.0  -t-  k 4-  8/c; 


(79) 


where  | SA: | « \k0  + K|,  k„  « K,  &>/ v.„  » k and  K — <o/v.„.  Substituting  Eq.  (79)  into 
(78)  we  find  that  the  spatial  growth  rate  is,  for  a highly  relativistic  beam, 
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The  result  (80)  has  been  previously  found  by  Sprangle  et  al.  (1975),  Kwan  et  al.  (1977), 
Kroll  and  McMullin  (1978),  and  Sprangle  and  Drobot  (1978).  This  result  is  found  under  the 
assumption  |8*|  < < 2k,  which  leads  to  the  constraint  upon  the  pump  amplitude 
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The  difference  between  the  phase  velocity  of  the  longitudinal  wave  and  the  axial  beam  velocity 


is 
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Using  this  result  in  the  efficiency  expression  of  Eq.  (46),  we  find  (Kwan  et  al , 1977;  Sprangle 
and  Drobot,  1978) 


(b)  Strong-Pump  Regime 

In  the  weak-pump  regime,  where  the  longitudinal  wave  satisfies  1 + x(<o,k)  = 0,  we  see 
from  Eq.  (74)  that  the  magnitude  of  the  scalar  potential  4>(z),  owing  to  charge  separation,  is 
much  greater  than  the  magnitude,  of  the  ponderomotive  potential  z ).  Thus  the  space-charge 

collective  effects  dominate  the  interaction.  When  the  pump  amplitude  is  sufficiently  high,  how- 
ever, the  ponderomotive  forces  may  modify  the  electrostatic  wave  dispersion  so  that 
Ix(«j.*)|  « 1.  In  this  case  the  space-charge  collective  effects  are  negligible.  Thus  to  consider 
the  strong-pump  regime  we  take 

|8Ar  | » 2k  (85) 

in  Eq.  (79).  Substituting  (79)  into  (78)  the  dispersion  relation  reduces  to 

eta  ! i.) 
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and  the  unstable  root  is  (Kroll  and  McMullin,  1978;  Sprangle  and  Drobot,  1978) 

hk~  ' (gftm >2/3*o-  (87) 

From  the  result  in  (87),  we  see  that  the  strong-pump  condition  (85)  is  the  inverse  of  (81),  i.e., 

0oi  » 0ch.  = 4U/y.?o),/J.  (88) 

The  phase  velocity  of  the  unstable  longitudinal  wave  is 
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From  Eqs  (46)  and  (48)  we  find  that  the  corresponding  efficiency  and  saturated  amplitude  are 
given  by 
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and 


|e|m,  = M£4yXi/4)l/3  , (91) 

respectively  (Sprangle  and  Drobot,  1978).  The  energy  spread  requirement  to  use  a cold  beam 
distribution  is  given  by 


(92) 


The  general  behavior  of  the  growth  rate  and  efficiency  as  a function  of  (pump  strength),  in 
both  the  strong  and  weak-pump  regimes,  are  shown  schematically  in  Figure  6.  In  the  weak 
pump  regime  the  spatial  growth  rate  is  proportional  to  the  pump  strength,  while  the  efficiency  is 
independent  of  the  pump.  For  the  strong  pump  case  both  the  growth  rate  and  the  efficiency 
vary  as  /3^3.  Higher  growth  rates  and  efficiencies,  therefore,  make  the  strong  pump  limit  a 
more  attractive  mechanism  for  generating  high  frequency  radiation.  In  Table  I we  summarize 
the  results  of  analyses  of  all  the  cases  1-3  for  a highly  relativistic  beam  (v.  = c). 


IV.  DISCUSSION:  A VISIBLE-LIGHT  FREE  ELECTRON  LASER 


In  this  section  we  apply  the  theory  to  a specific  example  of  a high-power  FEL  at  optical 
wavelengths.  We  first  note  two  considerations  important  for  high-current  beams:  (i)  the  effect 
of  the  beam  itself  on  the  pump  magnetic  field;  and  (ii)  the  induced  energy  spread  in  the  beam 
due  to  its  self  electrostatic  potential.  These  considerations  place  an  upper  limit  in  the  beam 
current. 

With  regard  to  (i)  we  note  that  the  beam  equilibrium  motion  in  the  external  pump  field 
contains  a transverse  velocity  component  0 ,dc.  The  beam  motion,  therefore,  induces  a spatially 
periodic,  diamagnetic  field  that  tends  to  oppose  the  applied  pump  field.  It  may  easily  be  shown 
that  with  an  equilibrium  distribution 

g„(a,fi,u)  = no~  8(a)8(/3)8(t/  — u0), 
the  diamagnetic  field  amplitude  B is  of  order 

I /I  ~ (2b0. 

where  £ is  the  beam  strength  parameter  defined  above:  £ = u>J{JT0ckJ.  Thus,  in  order  to 
neglect  the  perturbations  to  the  equilibrium  orbit  owing  to  the  diamagnetic  field,  we  require 
£ « 1. 

With  regard  to  point  (ii),  we  note  that  an  un-neutralized  electron  beam  has  an  inherent 
energy  spread  owing  to  the  variation  of  the  electrostatic  potential  within  the  beam.  For  a 
cylindrical  beam  of  radius  r„  the  potential  difference  between  the  central  axis  to  the  outer  edge 
of  the  beam  is 

- ir|<?  | n„r2. 

Associated  with  this  potential  difference  is  a kinetic  energy  change  across  the  beam,  given  by 

A£  = bym„c2  =■  |e|A</>. 

This  variation  in  kinetic  energy  may  be  interpreted  as  indicating  a thermal  energy  spread  in  the 
beam,  of  relative  magnitude 
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where  we  have  taken  y0  » 1.  For  sufficiently  high  beam  densities  this  inherent  energy  spread 

may  invalidate  the  cold-beam  approximation. 

As  noted  above,  the  strong-pump  scattering  require  (case  3b)  is  characterized  by  the 
highest  spatial  growth  rates  and  efficiencies,  and  is  thus  the  most  promising  case  for  application. 
Therefore,  we  shall  choose  parameters  appropriate  to  this  regime. 


We  take  a periodic  magnetic  pump  field  of  period  / * 0.75  cm  and  amplitude  B - 7.5  kG. 

The  electron  beam  current  /„  is  taken  to  be  2 kA,  with  beam  radius  0.1  cm  and  energy  45  Mev 
(y0  — 90).  The  transverse  velocity  induced  by  the  pump  is 

Vai  = Pai  c = \e\BJy „m0ck0  - 6.0  x I0'3c, 

where  k0  — 2n/l.  The  radian  plasma  frequency  of  the  beam  is  uih  = 2.0  x 10"  sec-1,  and  the 
pump  strength  parameter  f = Mt,/(y/yJ,ck0)  = 8.4  x 10-2.  The  longitudinal  Lorentz  factor 
y:o  — 80,  giving  p cr„  = 4(£/y?„)l/2  — 1.62  x 10-3  « verifying  that  these  parameters  are 
indeed  characteristic  of  the  strong-pump  regime.  The  wavelength  of  the  scattered  wave  is 
X,  — // 2y}„  = 0.58  /am,  in  the  visible  region  of  the  spectrum.  From  the  formulae  of  Case  3b, 
summarized  in  Table  1,  we  find  a spatial  growth  rate  |-ImX|  = JJ(.^pil±/4)2n  kn  = 3.6  x 10-2 
cm-1,  corresponding  to  an  e-folding  leng’K  of  28  cm.  Thus,  very  high  gains  may  be  achieved  in 
a single  pass  over  a total  distance  of  several  meters.  The  single-pass  efficiency  is  given  by 
V = (f/3ni/4)2/3  — 0.25%.  With  regard  to  the  constraints  on  the  current  imposed  by  the  neglect 
of  the  induced  diamagnetic  field  and  the  beam  quality  necessary  to  apply  the  cold-beam  approx-  , 

imation,  the  first  of  these  is  well  satisfied  as  £2  « 1,  and  the  second  is  marginally  satisfied  as 
A £/£„  = (£r0k0/ 2) 2 = 0.12%,  to  be  compared  with  the  efficiency  of  0.25%.  The  radiated 
power  £„  at  0.58  ^im,  is  given  by  /*,  = 17 /„£„  — 230  MW.  The  parameters  of  this  example  are 
summarized  in  Table  II. 

Finally,  we  note  that  it  is  in  principle  possible  to  increase  substantially  the  single-pass 
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decreasing  the  phase  velocity  of  the  pondermotive  wave  in  order  to  delay  trapping,  which 


saturates  the  process.  The  decrease  in  phase  velocity  can  be  achieved  by  adiabatically  decreas- 
ing the  wavelength  of  the  pump  field.  In  this  way  energy  will  be  extracted  from  the  beam  until 
a greater  degree  of  thermalization  is  attained,  with  correspondingly  higher  efficiency.  Extrapo- 
lating the  results  of  present  technology  to  the  FEL,  we  might  expect  the  combination  of  beam 
recovery  and  phase-velocity  tapering  to  yield  increases  of  roughly  an  order  of  magnitude  in  the 
single-pass  efficiency.  As  of  this  writing  (mid-1978),  however,  detailed  analyses  of  these  effects 
have  not  been  completed;  these  will  be  the  subjects  of  continued  investigations. 


V.  LABORATORY  EXPERIMENTS 


Data  on  laboratory  studies  of  stimulated  scattering  from  relativistic  electrom  beams  are 
summarized  in  Table  III.  The  five  experiments  listed  divide  naturally  into  two  regimes;  i)  the 
low-gain  regime  investigated  using  the  superconducting  linear  accelerator  at  Stanford  University 
(examples  1 and  2 in  the  table);  and  ii)  the  high-gain  collective  regime  investigated  using 
intense  relativistic  electron  beam  accelerators  at  the  Naval  Research  Laboratory  and  at  Colum- 
bia University  (examples  3,  4,  and  5 in  the  table). 

The  low-gain  regime  studies  were  characterized  by  relatively  high  electron  energies  (24-43 
MeV)  and  correspondingly  short  output  wavelength  (3-11  /im).  Line  broadening  of  spontane- 
ous emission  was  dominated  by  finite  length  effects  rather  than  by  spread  in  electron  energy 
since  !/L  > 2 EJE„  (see  Eq.  38).  The  current  was  relatively  small  (0.1-2. 6 A)  corresponding  to 
a low  value  of  amplitude  gain;  in  an  amplifier  experiment  (example  1),  G,  = 0.035  <<  1.  In 
an  oscillator  experiment  (example  2),  a peak  output  power  of  7 kW  was  realized  which 
corresponded  to  0.006%  of  the  energy  in  the  electron  beam;  average  output  power  was  0.36 
watts. 

In  contrast,  the  intense  beam  studies  had  lower  electron  energy  (.86-2  MeV)  and 
correspondingly  longer  output  wavelengths  (400-1500  /*m).  Line  broadening  was  dominated  by 
electron  energy  spread  since  2 E,J E„  > l/L.  (£,„  was  estimated  to  equal  the  electrostatic 
energy  shear  across  the  beam  ) The  current  was  much  larger  (5-40  kA)  making  possible  the  par- 
ticipation of  collective  electron  beam  modes  in  the  stimulated  scattering  process,  and  resulting 
in  large  gain;  in  one  superradiant  oscillator  experiment  (example  3),  the  measured  amplitude 
gain  was  G,  >2.  In  a second  superradiant  oscillator  experiment  (example  3),  the  strength  of 
the  pump  was  increased  to  ensure  saturation  of  wave  growth,  and  a peak  power  output  of  8 
MW  was  measured;  this  was  three  orders  of  magnitude  larger  than  the  saturated  peak  power  in 
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the  low-gain  oscillator  experiment,  corresponding  both  to  larger  peak  beam  power  and  to  a 
somewhat  larger  efficiency  (0.2%).  As  yet,  there  is  no  continuous  power  capability  in  the 
intense  beam  experiments,  the  experiments  being  characterized  by  a single  electron  pulse  of 
10-50  nsec  duration. 

In  the  intense  beam  experiments  two  types  of  pump  wave  were  investigated:  a magnetic 
wiggler  (examples  4 and  5)  as  was  used  in  the  low  gain  experiments  and,  alternatively,  a power- 
ful electromagnetic  wave  (example  3).  The  intense  beam  experiments  also  had  an  externally 
imposed  uniform  axial  magnetic  field  that  was  large  enough  so  that  magnetic  resonances  in  the 
output  power  were  observed. 

Some  appreciation  for  the  significance  of  the  experimental  results  listed  in  Table  III  may 
be  obtained  by  comparing  the  peak  powers  achieved  with  other  existing  capabilities  for  generat- 
ing high-power  coherent  radiation.  Two  distinct  technologies  are  currently  being  used. 
Conversion  of  electron  kinetic  energy  into  microwaves  by  coupling  beam  waves  to  electromag- 
netic waves  on  a slow-wave  structure  is  a highly  developed  technology  at  centimeter  and  mil- 
limeter wavelengths.  This  technology  is  characterized  by  numerous  device  types  which  can  be 
scaled  or  tuned  to  cover  a wide  frequency  range.  As  shown  by  the  solid  line  in  Fig.  7,  how- 
ever, the  available  power  decreases  rapidly  as  wavelength  is  shortened. 

At  optical  and  near  infrared  wavelengths,  the  technology  of  lasers  based  on  atomic  and 
molecular  transitions  has  been  successful  in  generating  very  high  power  using  a variety  of 
media.  The  number  of  atomic  or  molecular  transitions  which  result  in  laser  action  are  limited, 
however,  and  high  power  is  generally  available  only  in  narrow  bands.  Some  tunable  molecular 
systems  do  exist  (e.#.  dye  lasers)  but  these  are  restricted  in  their  power  output  because  large 
average  power  operation  results  in  decomposition  of  the  active  medium  (/>.,  the  dye). 
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In  contrast,  the  preliminary  results  obtained  with  the  free  electron  laser  point  to  the  possi- 
bility of  developing  high-power,  tunable  systems  operating  over  a wavelength  range  extending 
from  the  millimeter  to  the  optical.  The  data  points  indexed  as  3 and  4 in  Fig.  7 correspond  to 
the  intense  beam  scattering  results  listed  in  Table  III  as  examples  3 and  4,  respectively.  The 
low-gain  oscillator,  result  obtained  with  the  Stanford  superconducting  linear  accelerator  (exam- 
ple 2 in  Table  III)  is  shown  by  the  third  data  point.  A more  detailed  description  of  the  various 
stimulated  scattering  experiments  follows. 

a)  Experiments  in  the  low  gain  regime  at  infrared  wavelengths. 

The  experimental  arrangement  corresponding  to  example  1 in  Table  III  is  depicted  in  Fig. 
8.  The  output  from  a 50  kW,  10.6  jzm  laser  was  amplified  upon  interacting  with  a beam  of  24 
MeV  electrons  passing  through  a wiggler  magnetic  field  with  period  / - 3.2  cm  The  axial 
electron  energy  in  the  interaction  region  was  evidently  y.  = 39  to  satisfy  //2y.2, -A,;  this  may 
be  compared  with  a total  energy  of  y„  = 48. 

As  shown  by  Kroll  and  McMullin  (1978),  the  experimental  parameters  indicate  that  the 
interaction  took  place  in  the  low  gain,  tenuous  cold  beam  regime;  thus,  the  expressions  for 
Case  1 in  Table  1 are  applicable.  Gain  was  measured  to  be  a linear  function  of  current  for  5 < 
/ < 70  mA.  The  maximum  energy  gain  measurement  was  7%  corresponding  to  a stimulated 
emission  of  60  keV/electron  or  an  efficiency  of  0.25%.  The  measured  gain  is  somewhat  larger 
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in  case  1 of  Table  I when  this  is  multiplied  by  a filling  factor  of  2 x 10  2,  which  is  the  ratio  of 
the  electron  beam  area  to  the  area  of  the  applied  10.6  pm  radiation.  The  measured  efficiency 
approaches  the  theoretical  value  of  n/k„L  = 0.31%,  indicating  that  the  process  was  close  to 
saturation. 
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The  spontaneous  power  emitted  at  10.6  /am  is  plotted  as  a function  of  electron  energy  in 

I 

Fig.  9 together  with  gain  measured  when  the  C02  laser  radiation  was  applied.  The  full  width  at 
the  Me  points  of  the  spontaneous  power  curve  is  0.4%  which  agrees  well  with  a theoretical 
broadening  of  0.3%  due  to  the  finite  length  of  the  interaction  region.  The  shape  of  the  gain 
curve  showing  inverse  dependence  on  the  first  derivative  of  the  spontaneous  emission  curve 
also  agrees  with  theory. 

The  experimental  arrangement  corresponding  to  the  oscillator  in  example  2 of  Table  III  is 
shown  in  Fig.  10.  The  same  helical  magnet  undulator  was  used  as  with  the  amplifier  experi- 
ment shown  in  Fig.  8.  However,  energy  of  the  electrons  was  increased  to  43  MeV  so  that 
stimulated  emission  occurred  at  the  shorter  wavelength  of  3.4  /im.  The  electrons  are  formed 
into  bunches  only  1.3  mm  in  length  with  a spacing  between  bunches  of  approximately  25  m. 

The  spacing  between  the  resonator  mirrors  is  carefully  adjusted  so  that  the  round-trip  bounce 
time  for  the  3.4  /am  radiation  is  equal  to  the  time  interval  between  electron  pulses  entering  the 
cavity. 

The  emission  spectrum  of  the  laser  oscillator  above  threshold  is  shown  together  with  the 
spectrum  of  the  spontaneous  radiation  emitted  by  the  electron  beam  in  Fig.  11.  The  0.9% 
linewidth  of  the  spontaneous  emission  should  be  compared  with  theoretical  linewidth  of  0.6% 
due  to  the  finite  length  of  the  helical  undulator.  Lasing  in  the  optical  resonator  reduces 
linewidth  to  0.2%. 

The  peak  output  power  was  7 kW  corresponding  to  0.006%  of  the  beam  energy.  Assum- 
ing this  experiment  to  be  in  the  same  parameter  regime  as  the  amplifier  experiment  of  Elias  et 
al.  one  would  expect  that  the  applicable  theoretical  efficiency  would  again  be  n/k„L  - 0.3%. 

Perhaps  the  lower  experimental  efficiency  was  due  to  a choice  of  output  window  transmissivity 
that  was  not  optimum. 
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b)  Submillimeter  generation  in  the  high-gain,  collective  regime  with  an  electromagnetic  pump. 


The  experimental  arrangement  corresponding  to  example  3 in  Table  III  is  shown  schemat- 
ically in  Fig.  12.  A voltage  pulse  (-2  MV,  50  ns)  was  applied  to  a cold  cathode  immersed  in  a 
uniform  axial  magnetic  field  < B>  = 13kG.  A resultant  annular  electron  beam  was  injected 
into  a drift  tube  along  the  lines  of  <B>.  Electrostatic  forces  reduced  electron  kinetic  energy 
so  that  y„  - 3.9. 

At  one  point  in  the  drift  tube,  the  electron  beam  was  passed  through  a nonadiabatic  per- 
turbation in  the  applied  magnetic  field  which  converted  a large  fraction  of  the  electron  stream- 
ing energy  into  energy  transverse  to  the  axis  of  the  drift  tube.  The  electrons  with  large 
transverse  energy  then  reacted  unstably  with  an  electromagnetic  mode  of  the  drift  tube  (cyclo- 
tron maser  process  as  described  in  Hirshfield  and  Granatstein  (1977))  and  produced  a 2 cm 
pump  wave  having  a power  of  10-100  MW.  This  pump  wave  was  reflected  off  a metallic  plate 
and  travelled  back  up  the  drift  tube  toward  the  cathode.  When  the  pump  wave  encountered  the 
cold  streaming  electrons  near  the  cathode,  stimulated  scattering  occurred  resulting  in  backscat- 
tered  submillimeter  waves. 

The  submillimeter  output  radiation  was  monitored  through  an  array  of  metallic  tubes  pass- 
ing through  the  reflecting  plate.  Each  tube  passed  the  radiation  through  wave  filters  and  into  a 
pyroelectric  detector.  The  submillimeter  output  power  was  optimized  by  varying  both  the  axial 
position  of  the  magnetic  modulator  and  the  level  of  </#>.  The  total  submillimeter  radiation 
incident  on  the  metallic  reflector  was  shown  to  be  peaked  at  X,  = 400  and  have  a power  of 
1 MW. 

Dependence  of  the  submillimeter  radiation  on  the  position  of  the  magnetic  modulator  is 
shown  in  Fig.  13.  As  the  modulator  was  moved  downstream  from  the  cathode,  the  scattering 
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region  L2  was  lengthened  while  simultaneously  shortening  the  region  Z.,  in  which  the  pump 
wave  is  generated.  Thus,  one  expected  an  initial  increase  in  the  submillimeter  output  as  the 
scattering  became  stronger  followed  by  a decreasing  submillimeter  output  as  the  pump  wave  is 
weakened.  In  Fig.  13,  a plot  of  submillimeter  power  vs.  L,  shows  the  expected  behavior  with 
the  submillimeter  output  having  a maximum  at  L2  = 30  cm.  The  rising  portion  of  the  curve  in 
Fig.  13  has  a slope  indicating  an  amplitude  gain  for  the  stimulated  scattering  of  G,  > 2 or 
equivalently  a spatial  growth  rate  of  Im(A  ) > 0.02  cm  '1. 

Only  a lower  bound  is  determined  by  the  data  of  Fig  13  since  the  pump  wave  amplitude 
decreases  as  L2  is  made  larger.  This  lower  bound  on  spatial  growth  rate  may  be  compared  with 
the  theoretical  value  -ImU)  = 0.04  cm'1  calculated  from  the  expression  for  case  3a  in  Table  I, 
and  using  the  expression  for  fi„±  from  Sprangle,  et  at.,  (1975)  (viz.,  = 

k|6l>/-y,>m,)(&>„  - <n>/y„)  where  is  the  electric  field  and  < ft  > = |c|  <B>/m„c.  For 
the  parameters  of  the  experiment  in  Fig.  10.  we  calculate  0„±  = .02  and  = 0.4  while 
£■,;,/£„  = 0.04  and  £/y.„  - .13;  thus,  the  condition  /3,a  « /3,.r„  is  well  satisfied  while 

EJEo  « Zly-.n  is  marginally  satisfied.  Hence,  one  is  justified  in  using  the  theoretical  expres- 
sions for  the  high  gain  collective  scattering  regime  of  case  3a. 

Dependence  of  the  submillimeter  power  on  <B>,  the  solenoidal  magnetic  field  strength, 
is  shown  in  Fig.  14.  A maximum  is  seen  at  <B>  = 13  kG  with  a half-power  width  of  3 kG. 

The  pump  wave  power  at  Xs  = 2 cm  was  relatively  constant  in  the  range  10  kG  < <B>  < 16 

kG  so  that  the  behavior  depicted  in  Fig.  14  corresponds  to  resonant  interactions  in  the 
stimulated  scattering  process. 

The  experimental  configuration  in  Fig.  12  is  limited  by  the  fact  that  both  pump  wave  gen- 
eration and  stimulated  scattering  occur  in  the  same  electron  beam  This  does  not  allow  for 
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separate  optimization  of  each  process  and  complicates  interpretation  of  the  experimental  results. 
To  overcome  these  disadvantages,  the  new  experimental  configuration  in  Fig.  15  is  being 
assembled  at  NRL.  Two  synchronized  pulsed-power  systems  are  being  employed.  The  pulse- 
forming network  from  a 600  kV  flash  x-ray  system  (peak  power  - 3.6  GW)  will  be  used  to 
drive  an  S-band  magnetron  that  will  furnish  a multi-gigawatt  pump  wave  at  X„  = 9.4  cm. 

A larger  2.5  MV  accelerator  (peak  power  25  GW)  will  provide  the  annular  electron  beam 
(A,,  = 7 cm2)  in  which  stimulated  scattering  may  take  place.  The  backscattered  output  radiation 
is  expected  to  be  at  X,  = 1700  /xm.  Given  the  magnitude  of  microwave  power  available  from 
the  magnetron,  this  stimulated  scattering  experiment  is  expected  to  be  the  first  operating  in  the 
strong  pump  limit  (/3,a  > /3,„,).  The  expressions  for  case  3b  in  Table  I are  applicable  and 
using  them  one  calculates  a spatial  growth  rate  -1m (X)  = 0.4  cm'1  and  an  efficiency  rj  = 11%. 

The  investigations  of  stimulated  scattering  using  electromagnetic  pump  waves  are  of  prac- 
tical importance  since  pump  wavelengths  can  be  realized  which  are  much  shorter  than  those 
which  are  possible  with  a magnetic  wiggler.  Thus  for  a given  output  wavelength  there  is  a 
reduced  requirement  on  accelerator  voltage.  However,  as  is  apparent  from  the  discussion 
above,  the  production  of  strong  electromagnetic  pump  waves  is  far  more  difficult  technologi- 
cally than  the  production  of  a magnetic  wiggler  field. 

c)  Experiments  in  the  high-gain,  collective  regime  with  a magnetic  wiggler 

The  experimental  arrangement  corresponding  to  example  4 in  table  III  is  depicted  in  Fig. 
16  As  in  part  b)  above  an  intense  electron  beam  of  annular  cross-section  was  propagated  axi- 
ally down  a drift  tube  along  the  lines  of  a solenoidal  magnetic  field  <B>.  However,  in  this 
case,  a ripple  component  was  added  to  the  solenoidal  field  over  a distance  typically  L = 40  cm. 
The  radial  component  of  the  ripple  field  B„  could  be  varied  from  0 to  500  G,  and  ripple  periods 
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of  6 mm  and  8 mm  were  used.  From  the  parameters  in  table  III,  one  calculates  0(>I  = 0.01, 

“ 0-35,  and  £/y.„  - 0.03  while  EJE„  - 0.02.  Thus  fi„±  « /3ir„  is  well  satisfied,  but 
E/J E„ _<£/y;„  rather  than  E„J E„  « £/y;„.  One  therefore  expects  operation  to  have  been 
marginally  in  the  high-gain,  collective  regime. 

The  emission  spectrum  was  resolved  using  a grating  spectrometer  of  special  design.  A 
typical  spectrum  is  shown  in  Fig.  17.  As  shown  in  the  figure,  it  was  found  that  the  emission 
spectrum  consisted  of  peaks  at  two  distinct  wavelengths  which  were  interpreted  as  correspond- 
ing to  the  frequencies 

- (1  + P-Jy - otjy},17) 

and 

«S2“  (1  + P-Jy;o(vznk„  - |e|  <B>/mncyll). 

The  shorter  wavelength  peak  at  o>s,  would  then  correspond  to  a stimulated  scattering  process 
involving  a space  charge  beam  mode  as  an  idler.  The  longer  wavelength  peak  at  a»v2  would 
correspond  to  a scattering  process  with  a cyclotron  mode  idler.  Indeed  it  was  found  that  the 
longer  wavelength  peak  would  be  tuned  in  frequency  by  changing  <5>,  while  o,,  remained 
constant. 

Marshall  et  al.  (1977)  measured  power  in  the  wave  at  o>S2  to  be  8 megawatts.  This  was 
found  experimentally  to  represent  a saturated  power  level  (Gilgenbach,  1978)  and  corresponds 
to  an  efficieny  of  0.2%.  The  theoretical  analysis  presented  in  this  paper  does  not  apply  to  the 
case  of  a cyclotron  mode  idler.  The  power  measured  at  wv,  was  2 MW  but  this  was  not 
established  to  be  a saturated  value. 

The  experiment  of  Marshall  et  al.  (1977)  shown  in  Fig.  16  functioned  as  a superradiant 
oscillator.  Experiments  have  now  been  carried  out  (McDermott  et  al , 1978)  using  mirrors  to 
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form  a quasi-optical  cavity.  The  experimental  arrangement  corresponding  to  example  5 in 
Table  III  is  shown  in  Fig.  18.  The  separation  between  the  mirrors  was  1.5  m while  the  time 
duration  of  the  electron  beam  pulse  was  approximately  40  ns  so  that  radiation  would  execute 
only  4 round-trip  bounces  between  the  mirrors  while  the  beam  was  on.  Nevertheless,  evidence 
of  lasing  was  observed  with  coherence  of  the  radiation  measured  as  A A,/A,  = 2%  compared 
with  AA ,/A,  > 10%  for  the  superradiant  oscillator.  According  to  McDermott  et  al  (1978),  the 
stimulated  scattering  process  involved  a space  charge  wave  idler,  while  scattering  from  the 
cyclotron  idler  was  damped. 

VI.  CONCLUSIONS 

Free  electron  lasers  based  on  stimulated  scattering  from  intense  electron  beams  hold  the 
promise  of  making  available  very-high-power,  moderately  coherent,  widely  tunable  sources  over 
a wavelength  range  from  millimeters  down  to  optical  wavelengths  and  beyond.  This  paper  con- 
tains design  examples  of  a 1700  laser  with  peak  power  of  several  gigawatts  operating  at  11% 
efficiency  and  a 0.58  /im  laser  with  power  of  hundreds  of  megawatts  operating  with  0.25% 
efficiency. 

Clearly,  much  investigation  will  be  required  before  this  potential  can  be  fully  evaluated. 
The  experiments  which  have  been  carried  out  to  date  give  support  to  the  linear  theory,  but 
carefully  designed  experimental  investigations  of  the  nonlinear  regime  still  remain  to  be  exe- 
cuted. 

Also,  accelerator  technology  studies  are  required  to  evaluate  the  average  power  potential. 
In  connection  with  free  electron  lasers  based  on  linear  accelerators,  the  use  of  storage  rings  is 
under  study.  For  free  electron  lasers  based  on  intense  electron  beams,  there  is  great  interest  in 
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the  studies  which  have  been  initiated  on  repetitively  pulsing  intense  beam  accelerators  and 
recovering  energy  in  the  spent  electron  beam  through  depressed-collector  techniques. 

The  theoretical  studies,  while  presently  more  advanced  than  the  experiments,  are  also  far 
from  complete.  It  should  be  of  great  value  to  relax  some  of  the  model  idealizations,  and 
include  consideration  of  oblique  angle  scattering,  energy  shear  in  the  electron  beam,  the  pres- 
ence of  a uniform  magnetization,  and  finite  transverse  beam  geometry.  Most  important  of  all, 
the  study  of  axially  nonuniform  systems  (e.g.,  tapering  the  periodicity  of  the  magnetic  wiggler) 
may  indicate  efficiencies  greatly  enhanced  over  those  calculated  above  for  uniform  systems. 
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APPENDIX 


We  indicate  briefly  here  the  details  of  the  derivation  of  Eqs.  (24)  from  Eq.  (20),  within 
the  context  of  the  distribution  function  (23).  Substituting  the  expression  for  £(l>  given  in  Eq. 
06)  with  £(0)(a,  p,  u)  = n„ 8(a)  S(/3)  gn(u ) into  Eq.  (20)  and  integrating  by  parts  gives 
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In  order  to  evaluate  the  current  densities  in  (Al)  the  following  relations  are  needed: 


t a \ T ■ T ■ a \e\ 

p±(a,  P,  z ) = py  T ipr  -uT/jS — e ; 

c kn 


p(a,  P,  u,  z)  - u2  - a + Aoy(z)  - P + A „,.(z)  ; 

c c 


"01  “ “ f»„  HA.  “ M BJm0c ; 

C Kn 


P±(a,  P,  u,  z)  J ^ Wqi  e*’1'"  ', 

P=(a.  P.  u,  z)  | - y/u2  - u2±  - u.(u)  - -y(u)  m„  v.(u); 

a— p— 0 

|ir  + iw]  '* <o • *• !>  .-L  * lol; 

-U  ‘ l“) 


48 


d 3 
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With  the  above  relations  it  becomes  straightforward  to  derive  Eqs.  (24). 
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Table  I — Summary  of  FEL  Regimes 


mnc  , E/if  * y :0 y 0 fttn c V,/n 


Table  II 

Illustration  of  Visible  Radiation  Source  (X,  — 0.58  ^m) 
(Strong  Pump  Scattering  Regime) 

Magnetic  Pump  Field  Parameters 


Pump  Wavelength:  / — 0.75cm 


Pump  Magnetic  Field  Amplitude: 

B„  = 7.5  KG 

Electron  Beam  Parameters 

£-Beam  Current: 

/„  - 2k A 

c'-Beam  Energy  (y„  = 90): 

£„  - 45  MeV 

£-Beam  Power: 

P„  = 90  GW 

£-Beam  Radius: 

r„  — 0.1  cm 

Axial  Beam  Gamma: 

H 

oo 

o 

Energy  Spread  Due  to  Self 

Beam  Potential: 

A£/£„  = 0.12% 

Scattered  Visible  Radiation  Parameters 

Radiation  Wavelength: 

X,  = l/(2y}„)  - 

e-Folding  Length  of  Radiation: 

£,.  = 28  cm 

Single  Pass  Efficiency: 

V - 0.25% 

Radiation  Power: 

Pt  “ 230  MW 

0.58  /im 
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WIGGLING 


PERIODIC  CIRCULARLY  POLARIZED 
MAGNETIC  PUMP 


-INTERACTION  REGION  - 

Figure  1 — Schematic  of  the  free-electron  laser  model.  The  unmodulated  electron  beam  enters 
the  interaction  region  from  the  left.  The  pump  field  builds  up  adiabatically  and  reaches  a con- 
stant amplitude  for  r > 0. 
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(y/y.=1.4) 


F W0.  k0L,/?1yz) 


k0L=2rr 

(1  PERIOD) 
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Figure  2 (Continued)  — The  gain  function  F(90,k0L,  p±  yz)  corresponding  to  the  low-gain, 
tenuous  cold  beam  scattering  limit.  The  interaction  length  for  (a)  is  160  pump  periods  and  for 
(b)  1 pump  period.  Note  that  in  (b)  the  slight  asymmetry  with  respect  to  90  is  due  to  the  con- 
tributions of  the  terms  sin2#,,/#,,  and  sin2(#„  + k„L/2)/(90  + k0L/2). 
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PHASE  VEL.  OF  LONGITUDINAL 
WAVES 


Figure  3 — Pictorial  representation  of  phase  velocities  of  the  longitudinal  waves  and  electron 
beam  distribution  function.  When  the  phase  velocities  lie  far  outside  the  beam  distribution,  the 
beam  may  be  considered  cold. 


Figure  4 — Schematic  of  quasilinear  plateau  formation  in  saturation  of  relativistic  stimulated 
Compton  scattering  limit. 
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Figure  5 — Coupled  dispersion  diagrams  of  the  beam  waves  and  transverse  electromagnetic 
waves.  Instability  occurs  near  the  intersection  between  the  slow  beam  made  and  transverse 
electromagnetic  made. 
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Figure  6 — Schematic  of  spatial  growth  rate  and  efficiency  in  the  high-gain,  dense 
cold  beam  limit  as  a function  of  the  pump  field  strength  characterized  by  pol. 
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Figure  7 — A Comparison  of  Technologies  Available  for  the  Generation  of  High-Power 
Coherent  Radiation.  The  three  data  points  correspond  to  the  preliminary  coherent  scattering 
results  and  are  numbered  as  in  Table  III. 
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Figure  8 — Experimental  arrangement  in  free  electron  amplifier  study  in  the  low-gain  regime. 
Fig.  8 is  reproduced  from  Elias  et  al  (1976). 
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Figure  9 — (a)  Spontaneous  power  at  10.6  fim  as  a function  of  electron  energy,  and  (b)  the 
gain  in  energy  of  10.6  Mm  radiation  imposed  from  the  C02  laser  as  a function  of  electron 
energy.  Fig.  9 is  reproduced  from  Elias  et  al  (1976). 
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Figure  10  — Experimental  arrangement  in  free  electron  laser  oscillator  study  in 
the  low-gain  regime.  Fig.  10  is  reproduced  from  Deacon  et  al  (1977). 
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ABOVE 

THRESHOLD 


(a) 


Figure  1 1 — Emission  spectra  of  laser  oscil- 
lator (a)  above  lasing  threshold  (b)  spon- 
taneous radiation  emitted  by  electron  beam. 
Fig.  11  is  reproduced  from  Deacon  et  al., 
(1977). 


0.7  my 


INSTRUMENT 

WIDTH 


■ 8 mu 


-i_L 


(bl 


3.417  M 


REGION  GENERATION 

REGION 


Figure  12  — Experimental  arrangement  in  study  of  stimulated  scattering  in  the  high-gain,  col- 
lective regime  with  electromagnetic  pump  wave  (Granatstein  et  al,  1977)  The  VEBA  electron 
accelerator  system  is  described  in  Parker  and  Ury  (1975). 
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Figure  14  — Submillimeter  output  radiation  at  400  / 
field  (Granatstein,  et  al,  1977). 
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Figure  15  — Planned  Stimulated  Scattering  Experiment  with  Electromagnetic  Pump  Wave. 
Operation  in  the  high-gain,  strong-pump  regime  is  expected  (Case  3b  in  Table  I). 


SOLENOID 


SPECTROMETER 

Figure  16  — F.xperimental  arrangement  in  study  of  stimulated  scattering  in  the  high-gain,  col- 
lective regime  with  magnetic  undulator  pump.  Fig.  16  is  reproduced  from  Gilgenbach  et  al 
(1978). 
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Figure  17  — Spectrum  of  output  radiation.  Each  peak  arbitrarily  normalized  to 
100  on  vertical  scale.  / = 8 mm,  y20  = 1.8,  <B>  - 7.6  kG  for  peak  at  tus2, 
n0  = 3 x 10"cm-3.  Fig.  17  is  reproduced  from  Gilgenbach  et  al.,  (1978). 
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Figure  18  — Experimental  arrangement  in  study  of  high-gain  free  electron  laser  with  quasi- 
optical  cavity.  A/,  and  M2  designate  mirrors.  Fig.  18  is  reproduced  from  McDermott  et  al 
(1978). 
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